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Volume  I of  of  this  rep<N’t  presents  the  results  of  a research  program  to  evaluate  and  develop  water-wave  theories 
for  engineering  application.  Volume  11  presents  wave  tables  developed  for  preliminarv-  design  in  offshore  problems. 


^ Volume  I deschlies:  (l)an  evaluation  of  the  degree  to  which  various  available  wave  theories  satisfy  the 
nonlinear  water-wave  mathematical  formulation  and  (2)  a comparison  of  water  particle  velocities  measured  in  the 
Uborattwy  with  those  predicted  by  ■ number  of  available  wove  theories.  The  results  indicated  that  Dean  s 
.S'tream-Fum’tion  Wave  Theory  (prcjvided  generally  better  agreement  with  bf>lh  the  mathematical  formulation  and 

((^tinued) 
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the  laboratory  data.  Volume  I abo  includes  a number  of  examples  illustrating;  the  applieatif>n  of  the  wave  tables 
(described  below)  to  offshore  de^gn  problems.  * — ^ 


•Xy  Baaed  on  the  evaluation  phase  described  above,  a set  of  wave  tables  was  deveio|>ed  and  is  pn^sented  as  Volume  II. 
The  tables  consist  of  dimensionless  quantities  which  describe  the  kinematic  and  dynamic  fields  of  a 
two^limensional  progressive  water  wave.  In  addition,  quantities  are  included  which  are  directly  applicable  to 
frequently  required  design  calculations  and  also  parameters  which  should  be  of  interest  to  the  researcher  and 
scientist.  


PREFACE 


Thr  U.S.  Army  Coastal  Engim-erinp  Rest-arch  Center  (CERC)  is  puhlisliing  this  report  and 
eolh-ction  of  Siream-fnnetion  tables  to  assist  coastal  enpneers  in  the  design  of  structures 
vulnerable  to  wave  forces.  The  report  was  prepared  by  R.  G.  Dean,  College  of  Engineering, 
University  of  Florida. 

The  author  acknowledges  the  help  of  other  peo|)le  in  the  preparation  of  this  report. 
E.  Olsen  aiio  B.  Beechley  made  many  of  the  computer  calculations  and  organized  the 
results.  R.  A.  Dalyrymple  assisted  in  the  later  pha.ses  by  applying  some  of  the 
.Stream-function  computer  techniques  which  he  has  developed.  M.  P.  O’Brien,  senior  civilian 
member  of  the  Coastal  Engineering  Research  Board  offered  constructive  suggestions  and 
dist  ussions.  J.  Richard  Weggel,  Special  Assistant  to  the  Commander  and  Director,  CERC, 
made  several  useful  suggestions  which  simplified  the  application  of  the  tables  to  design 
problems.  D.  Lee  Harris,  the  contract  monitor  for  CERC,  provided  encouragement  and 
useful  comments  regarding  the  format  of  the  tables.  Reviews  of  an  early  draft  by  persons 
actively  involved  in  offshore  design  or  wave  theory  studies  assisted  in  preparation  of  the 
n-j)ort.  The  constructive  comments  and  suggestions  from  these  reviews  are  greatly 
appreciated. 

This  report  was  prepared  for  the  Coastal  Engineering  Research  Center  under  Contract 
DACW  72-67-C-0()09  with  the  University  of  Florida. 

NOTE:  Comments  on  this  publication  are  invited. 

Approved  for  publication  in  accordance  with  Public  Law  166,  79*^  Congress  approved  31 
July  1949,  as  supplemented  by  Public  Law  172, 88*^  Congress,  approved  7 November  1963. 
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Symbol  Description 

BBC  Itottom  boundary  condition,  defined  by  Equation  (10) 

C Wave  celerity 

Cj)  coefficient 

Group  velocity 

Cj;  Dimensionless  group  velocity,  defined  by  Equation  (42) 

C^f  Inertia  coefficient 

D Pile  diameter 

DFSBC  Dynamic  free  surface  boundary  condition,  defined  in  Equation  (12) 

DFSBP  Dynamic  free  surface  breaking  parameter  defined  in  Equation  (49) 

~D  Subscript  “D”  denoting  “design”  value,  also  drag  component  of  force 

and  moment 

Cj  Dynamic  free  surface  boundary  condition  error,  utibzed  by 

Cbappelear  and  Von  Schwind  and  Reid  defined  in  Equation  (17a) 

E|  Root  mean  square  error  in  kinematic  free  surface  boundaiy^  condition 

Ej  Root  mean  square  error  in  dynamic  free  surface  boundary  condition 

E Mean  square  error  in  dynamic  free  surface  boundary  condition 

Fp  Drag  force  component 

F'd  Dimensionless  drag  force  component,  defined  by  Equation  (25) 

Fj  Inertia  force  component 

F'l  Dimensionless  inertia  force  component,  defined  by  Equation  (26) 

Ffg  Total  energy  flux  in  direction  of  wave  propagation,  per  unit  width 
Dimensionless  form  of  Fjp,  defined  by  Equation  (41) 
f Momentum  flux  in  direction  of  wave  propagation 

Fj„^  Dimensiordess  form  of  F^^,  defined  by  Equation  (44) 

F^^  Momentum  flux  transverse  to  direction  of  wave  propagation 
FJ„^  Dimensionless  form  of  F^^,  defined  by  Equation  (45) 
g Gravitational  constant 

h Water  depth 

h'  Freeboard  used  in  establishing  deck  elevation 

h^^  Breaking  water  depth 
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Symbol  Description 

H Wave  height 

Breaking  wave  height 
Deepwater  wave  height 
j Index  used  in  summation 

J Maximum  value  of  j in  summation 

KE  Kinetic  energy  of  waves 

KE*  Dimensionless  form  of  KE,  defined  by  Equation  (39) 

KFSBC  Kinematic  free  surface  boundary  condition,  defined  in  Equation  (11) 
KFSBP  Kinematic  free  surface  breaking  parameter,  defined  by  Equation  (48) 
L Wavelength 

L'  Dimensionless  form  of  L,  defined  by  Equation  (37) 

Lg  Small-amplitude  deepwater  wave  length  = gT^/(2ir) 

~L  Subscript  “L”  denoting  “lower” 

M^)  Drag  moment  component 

Mj)  Dimensionless  form  of  , defined  in  Equation  (27) 

Inertia  moment  component 

M|  Dimensionless  form  of  Mj,  defined  in  Equation  (28) 

n Index  used  in  summation 

NN  Order  of  wave  theory 

p Pressure 

p^  Atmospheric  pressure 

p^  Dynamic  component  of  wave  pressure 

p^  Dimensionless  form  of  p^,  defined  in  Equation  (29) 

Q Bernoulli  term,  defined  in  Equation  (8) 

Q Average  value  of  Q 

S Vertical  coordinate,  referenced  to  bottom,  positive  upwards 

t Time  coordinate 

T Wave  period 


xii 


SYMBOLS— Continued 


Symbol  Description 

u Horizontal  component  of  water  particle  velocity 

o'  Dimensionless  form  of  u,  defined  in  Equation  (21) 

0^^  Measured  horizontal  component  of  water  particle  velocity 
uj-  Theoretical  horizontal  component  of  water  particle  velodt'' 

Subscript  “U”  denoting  upper  value 
w Vertical  component  of  water  particle  velocity 

w'  Dimensionless  form  of  w,  defined  in  Equation  (22) 

^L,U  Weighting  coefficients,  defined  by  Equation  (50) 

X Horizontal  coordiiute 

X Stream -function  coefficients 

* Vertical  coordinate,  referenced  to  still  water  level,  positive  upwards 

a Wave  crest  alignment  relative  to  bottom  contours 

Oq  Deepwater  value  of 

Wave  crest  alifftment  at  breaking  conditions 
0,  Kinematic  free  surface  breaking  parameter,  defined  by  Equation  (48) 

Pj  Dyiumic  free  surface  breaking  parameter,  defined  by  Equation  (49) 

7 Specific  weight  of  water 

Cj  Distributed  error  in  kinematic  free  surface  boundary  condition 

fj  Distributed  error  in  dynamic  tree  surface  boundary  condition 
e'j  Dimensionless  form  of  e, , defined  by  Equation  (36) 

V Water  surface  displacement 

T}'  Dimensionless  form  of  t),  defined  by  Equation  (30) 

6 Phase  an^e 

s Numerical  corutant,  3.14159  .... 

p Maas  density  of  water 

o Standard  deviation 

^ Velocity  potential 

^ Stream  function 

Stream-function  value  evaluated  on  free  surface  of  an  “arrested"  wave 
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EVALUATION  AND  DEVELOPMENT  OF  WATER  WAVE  THEORIES  FOR 
ENGINEERING  APPLICATION 


by 

R.  G.  Dean 
1.  INTRODUCTION 

The  following  were  the  primary  goals  of  the  research  reported:  (l)for  given  wave 
conditions,  to  establish  a rational  basis  for  selection  of  one  of  the  numerous  available 
progressive-water-wave  theories  and  (2)  to  tabulate  the  most  appropriate  wave  theory  or 
theories  in  a form  convenient  for  preliminary  design  use.  The  main  emphasis  has  been  an 
attempt  to  assist  the  engineer  in  his  selection  and  a})plication  of  wave  theories  in  marine 
design  j)roblems.  The  rcsi’arch  has  proceeded  in  several  distinct  phases  which  are  described 
briefly  Indow. 

An  early  phase  of  the  research  was  related  to  evaluating  the  analytical  validity  of 
water-wave  theories;  that  is,  the  degree  to  which  the  various  available  theories  satisfy  the 
etpiatioris  constituting  the  mathematical  formulation.  The  results  of  this  phase,  first 
published  in  September,  1968  (Dean,  1968a),  established  that,  of  the  eight  theories  included 
in  the  study,  the  Stream-function  fifth-order  provided  the  best  fit  over  a wide  range  of  wave 
conditions.  For  very  shallow  water  waves,  the  Airy  and  first-order  Cnoidal  theories  provided 
the  best  fit.  However,  because  the  Stream-function  theory  can  be  extended  to  quite  high 
orders,  it  was  expected  that  it  would  provide  the  best  fit,  even  for  most  shallow  water  wave 
conditions.  Based  on  the  results  of  this  phase,  the  following  phases  concentrated  on  further 
exploration  and  development  of  the  Stream-function  theory  for  engineering  application. 

The  second  phase  represented  an  examination  of  near-breaking  wave  conditions  using  the 
.Stream-function  theory  (Dean,  1968b).  This  problem  is  complicated  because  breaking 
conditions  represent  a mathematical  as  well  as  a hydrodynamic  instability,  and  therefore  the 
computational  aspects  are  not  straightforward.  The  results  of  this  study  indicated  that  of 
the  two  stability  criteria,  the  kinematic  criterion  rather  than  the  dynamic  criterion  governs 
at  breaking.  It  was  also  found  that  near  breaking  the  pressure  distribution  was  hydrostatic 
rather  than  characterized  by  a zero  pressure  gradient  as  predicted  by  .some  other  studies. 
The  complexities  of  the  numerical  computations  led  to  an  attempt  to  establish  the  breaking 
index  for  only  three  relative  water  depths  (shallow,  transitional  zone  and  deep).  It  was 
found  that  for  shallow  and  deepwater  waves,  the  breaking  heights  established  from  the 
.Stream-funcfi«>n  wave  theory  were  up  to  28  percent  higher  than  those  established  earlier  by 
other  investigations.  For  transitional  depth  conditions,  however,  the  breaking  heights 
determined  in  the  study  agreed  well  with  those  of  earlier  investigations. 

The  third  phase  of  the  investigation  (Dean  and  LeMehaute,  1970)  was  related  to  the 
“experimental  validity  of  water  wave  theories”  as  compared  to  “analytical  validity.”  The 
motivation  of  this  phase  was  the  recent  publication  of  a fairly  comprehensive  set  of 
measurements  of  water  particle  velocities  for  shallow  water  waves  and  comparison  witi;  a 
number  of  wave  theories  by  LeMehaute,  Divoky,  and  Lin  (1968);  a comparison  with  the 


Stream-function  theory  was  therefore  conducted  as  a part  of  tlie  pr«;sent  study.  On  an 
overall  basis,  tlie  Stream-function  theory  provided  a significantly  better  fit  to  the  measured 
water  particle  velocities  than  the  other  theories.  The  standard  deviation  between  the 
measured  and  Stream-function  representations  was  0.17  foot/second  as  compared  to 
0.24  foot/second  for  the  theory  providing  the  next  best  fit.  The  primary  significance  of  this 
phase  of  the  study  is  that  the  wave  conditions  are  in  the  .shallow-water  region  where  theories 
other  than  the  Stream-function  would  be  exjtected  to  provide  better  comparisons  with 
measurements.  Although  this  favorable  comparison  Ls  not  taken  as  demonstration  of  the 
superiority  of  tlie  Stream-function  for  all  wave  conditions,  the  results  were  very  encouraging 
and,  to  some  extent,  surprising. 

TTie  final  phase  of  the  investigation  has  been  the  development  of  a computer  program  to 
tabulate  wave  quantities  that  would  be  of  value  to  engineers  in  design,  and  that  would  also 
be  valuable  to  |)ers(yis  concerned  with  the  further  development  and  improvement  of 
water-wave  theories.  During  the  development  of  the  tables,  it  has  been  found  that  more 
meaningful  information  than  originally  anticipated  could  be  presented. 

In  the  early  phases  of  tliis  study,  dimensional  variables  [i.e.,  water  depth/(wave  period) 
and  wave  height/(wave  period)]  were  used  to  characterize  the  wave  conditions  (Dean, 
1968b).  This  feature  will  be  evident  in  the  description  of  some  of  the  results.  In  the  latter 
phases  of  the  study,  it  was  decided  to  characterize  the  wave  conditions  by  the  following 
dimensionless  quantities;  h/L^,  and  H/L^^,  where  h,  H and  represent  the  water  depth, 
wave  height,  and  small-amplitude  deepwater  wavelength,  respectively.  The  tables  are 
develof)ed  for  40  cases  of  (h/L^ , H/L^ ). 

The  results  of  the  research  are  presented  in  two  volumes.  The  present  report  (Volume  I) 
documents  the  research  results  and  describes  the  wave  tables  and  their  application. 
Volume  II  presents  the  wave  tables  that  have  been  developed  for  40  cases  encompassing 
most  conditions  encountered  in  engineering  design. 

It  .should  be  noted  that  all  of  the  available  wave  theories  have  not  been  included  in  the 
comparisons  described  earlier.  Some  of  the  theories  omitted  were  developed  during  the 
period  of  this  research;  some  have  been  available,  but  were  not  compared,  usually  becau.se 
they  are  not  employed  extensively  for  engineering  purposes. 

II.  STREAM  FUNCTION  WAVE  THEORY 

Introduction 

At  an  early  stagt;  of  the  research,  the  study  indicated  that  the  Stream -function  theory 
generally  provided  a better  fit  to  tlie  boundary  conditions  and  also  to  available  laboratory 
measurements.  The  study  therefore  develo|)ed  into  an  effort  to  explore  and  develop  the 
Stream-function  theory  for  engineering  application.  Before  presenting  this  work,  the  basis 
for  the  Stream-function  theory  will  be  described  in  some  detail  in  an  attempt  to  define  the 
similarities  with  and  differences  from  other  theories.  It  should  be  noted  that  there  are  two 
representations  of  the  Stream-function  theory:  (l)for  a given  wave  height,  H,  water 
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depth,  h,  and  wave  period,  T,  a (symmetrical)  representation  can  be  developed  to  describe 
the  kineniatii^i  and  dynamics  of  the  motion  and  (2)  for  a given  measured  water  surface 
displacement,  Tj(t)  representing  a single  oscillation  (e.g.,  trough-to-trough),  a representation 
can  be  determined  which  completely  defines  the  kinematics  and  dynamics  of  the  wave 
motion.  The  first  case  is,  of  course,  of  more  interest  to  designers;  in  another  application,  the 
second  case  has  been  employed  for  the  analysis  of  hurricane-generated  wave  and  wave-force 
data.  Only  the  first  mode  has  been  explored  under  the  present  study. 

Formulation 

The  water-wave  phenomenon  of  interest  here  can  be  idealized  as  a two-dimensional 
boundary  value  problem  of  ideal  flow.  The  a.ssumption  of  ideal  flow  is  essential  to  a 
mathematical  formulation  that  can  be  readily  solved  by  known  techniques.  Figure  1 defines 
terms  employed  in  the  formulation. 


Figure  1 . Definition  sketch,  progressive  wave  system 

Differential  Equation 

Ideal  flow  incorporates  the  assumptions  of  an  incompressible  fluid  and  irrotational 
motion.  For  pressures  normally  exjierienced  in  progressive  water-wave  motions,  the 
incompressibility  assumption  can  be  shown  to  be  valid.  Shock  pressures  due  to  a wave 
breaking  against  a seawall  may  be  an  important  exception;  however  these  are  not 
encompassed  by  the  results  of  this  research.  The  assumption  of  irrotational  flow  may  be 
questioned.  IVobably  the  best  reason  for  this  assumption,  at  this  stage,  is  that  it  allows 
formulation  of  a boundary-value  problem  that  can  be  solved  in  an  approximate  manner.  The 
solutions  can  then  be  compared  with  measurements  to  determine  the  apparent  need  for  the 
refinement  to  include  a nonzero  rotation. 

The  differential  equation  (DE)  for  two-dimensional  ideal  flow,  the  Laplace  equation,  can 
be  presented  in  terms  of  either  the  velocity  potential,  0 or  stream  function,  0, 

7^0  = 0 
* 0 
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(1) 

(2) 


where,  in  two  dimensions 


7=^ 


-3'  ^9* 

= * TIz 


(3) 


and  0 and  0 are  defined  in  terms  of  the  velocity  components  u and  w (see  Figure  1)  as: 


u 


w 


(4) 


Boundary  Conditionn 

Two  types  of  boundary  conditions  are  required  on  the  upper  and  lower  surfaces  (for  the 
present  study,  it  will  be  assumed  that  the  depth  is  uniform).  The  kinematic  boundary’ 
condition  applies  to  both  boundaries,  and  simply  requires  that  the  components  of  flow  at 
these  boundaries  be  in  accordance  with  the  geometry  and  motion  (if  any)  of  the  boundaries. 
This  condition  can  be  stated  as  follows: 

Bottom  Boundary  Condition  (BBC) 

w = 0,  z = -h  (5) 

Kinematic  Free  Surface  Boundary  Condition  (KFSBC) 

It  '^Ix  = W'  z = n(x,t)  (6) 

Dynamic  Free  Surface  Boundary  Condition  (DFSBC) 

The  remaining  free  surface  boundary  condition,  the  so-called  dynamic  free  surface 
boundary  condition  (DFSBC),  requires  that  the  pressure  immediately  below  the  free  surface 
be  uniform  and  equal  to  the  atmospheric  pressure,  p . 

■g'^"  constant  5 q',  z = n(x,t)  (7) 

In  the  above  formulation,  it  is  tacitly  assumed  that  surface  tension  effects  are  negligible.  It  is 
customary  to  incorporate  the  atmospheric  pressure  term  into  the  constant,  Q',  to  yield  a 
new  constant,  Q 

n + ^ (u^  + w*)  - i|i  = Q (8) 
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Ill  thr  furnuilatioii  prfsi'nled.  no  nMjuin'ments  liavc  brcti  |)lar(;d  on  the  permanence  of 
uvji’e  form:  that  is,  the  wave  eonld  ehanpe  form  as  it  propagates  due  to  the  relative  motion 
and  interference  of  components  propagating  with  various  phasi-  sjiecds.  The  treatment  of 
this  gi'iieral  problem  including  llie  nonlincaritics  is  complex,  and  was  not  the  subject  of  this 
research.  Rather,  in  the  presr-nt  investigation,  it  is  assumed  that  the  wave  propagates  witli 
constant  sjKcd,  C,  and  without  change  of  form.  It  i.s  then  possible  to  choo.se  a coordinate 
system  propagating  with  the  speed  of  and  in  the  same  direction  as  the  wave,  and  relative  to 
this  coordinate  system  the  motion  does  not  change,  and  is  therefore  steady.  The  time 
de|)«-ndencv  in  the  formulation  vanishes,  th<’  horizontal  velocity  component  with  rcsjiect  to 
the  moving  coordinate  system  is  u-(l;  and  the  formulation  may  be  .summarized  as: 


l)E:  V^4)  = = 0 


(9) 


Boundary 

Conditions 


BBC:  w = 0,  z = -h  (10) 

kfsfk::  ^ = n (X)  (ii) 

DF.'iBC:  n + ^ ((u  - C)^  + W^)-  = Q,  z = n(x)  (12) 

Motion  is  periodic  in  x with  spatial  jieriodicity  of  the  wavelength,  L.  (13) 


To  avoid  misimpressions  about  the  assumptions  and  formulation  presented  here  and  those 
employed  in  other  investigations  of  nonlear  waves,  it  is  noted  that  the  formulation 
incorporating  the  assumption  of  projiagation  without  change  of  form  is  common  to  the 
development  of  all  the  following  nonlinear  water  wave  theories: 

Stokes  2"*^,  and  higher  order  wave  theories 

t.noidal  l*t  and  2"*^  order  theories  by  e.g.,  Keulcgan  and  Patterson  (1940), 
and  Laitone  (1960) 

Sobtary  wave  theory,  1**  order  by  Boussinesq  (Munk,  1949) 

Solitary  wave  theory,  2^*^  order  by  McCowan  (Munk,  1949) 

Stream-function  wave  theory  by  Von  Schwind  and  Reid  (1972) 

To  reiterate,  analytical  mlidily  will  be  based  on  the  degree  to  which  a theory  satisfies  the 
boundary -value  problem  formulation.  Equations  (9)  through  (13).  If  a theory  could  be 
found  that  provided  exact  agreement  to  the  formulations,  then  the  analytical  validity  would 
be  jierfect.  There  is  no  guarantee  that  good  analytical  validity  ensures  that  a theory  will 
provide  a gtMid  representation  of  the  natural  phenomenon,  because  implicit  in  the 
formulation  are  the  a.ssumptions  that  capillary  and  rotation  forces  and  other  effects  are 
negligible.  Experimental  validity  will  be  based  on  the  agreement  between  wave  theories  and 
measured  data. 
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The  Stream-function  Solution 

For  the  formulation  expressed  in  Eciuations  (9)  through  (13),  a Stream-function  solution 
may  be  expressed  as: 

i)>  {x,2)  = ^ z + ' ^ X{n)  sinh  (h  + z)  cos  xj  (14) 

Evaluating  this  expression  on  the  free  surface,  i.e.,  setting  z = tj,  we  find 

n = ^ 'I'r,  - ^ X(n)  sinh  (h  -t-  n)j  cos  ^ x (15) 

where  NN  represents  the  order  of  the  representation,  i.e.,  the  number  of  terms 
contributing  to  the  series  expression,  represents  the  (constant)  value  of  the 
Stream-function  on  the  free  surface,  L is  the  (undetermined)  wavelength,  and 
the  X(n)  represent,  at  this  stage,  undetermined  coefficients. 

For  particular  wave  conditions,  it  is  regarded  that  the  wave  height,  period,  and  water 
depth  are  specified.  Equation  (14)  exactly  satisfies  the  governing  differential  equation  and 
the  bottom  and  free  surface  kinematic  boundary  conditions  for  arbitrary  values 
of  L,  and  the  X(n)  coefficients.  The  Stream-function  expression  is  also  periodic 
in  X with  wavelength,  L.  The  only  remaining  boundary  condition  is  the  dynamic 
free-surface  boundary  condition;  the  parameters  L and  the  X(n)’s  are  to  be  chosen  such 
that  this  boundary  condition  is  best  satisfied  for  a specified  wave  height. 

The  procedure  for  determining  the  unknown  parameters,  which  can  be  considered  as  a 
nonlinear  numerical  perturbation  procedure,  is  presented  in  Appendix  I. 

in.  EVALUATION  OF  VALIDITIES  OF  WAVE  THEORIES 

Introduction 

As  discussed  earlier,  there  arc  two  types  of  validity  that  were  examined.  “Analytical 
validity”  is  based  on  the  degree  to  which  a theory  satisfies  the  governing  equations  (of  the 
boundary  value  problem  formulation).  Good  analytical  validity,  however,  does  not 
necessarily  imply  good  representation  of  the  natural  phenomenon.  Experimental  validity  is 
based  on  the  agreement  between  a theory  and  measurements.  To  date,  some  reasonably 
good  laboratory  data  are  available,  and  at  least  two  field  measurements  of  water  particle 
velocities  are  reportedly  underway  (as  of  1972)  in  the  petroleum  industry,  and  hopefully, 
will  be  available  within  the  next  few  years. 

Diacuadon  of  Differences  Between  Stream-function  and  Other  Wave  Theories 

Later  in  this  section,  it  will  be  shown  that  the  Stream-function  theory  provides  a better 
fit  than  other  theories  to  the  boundary  conditions  and  also  provides  a better  fit  to 
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labiiratory  measurements  of  water  particle  velocities.  It  is  therefore  worthwhile  to  compare 
some  of  the  inherent  features  of  the  Stream-function  and  other  theories.  Although  it  is 
difficult  to  discuss  all  other  theories  in  general  statements,  an  attempt  will  be  made  to 
present  the  more  .significant  representative  differences. 

Consider,  as  an  example,  the  Stokes  higher  order  wave  theories.  The  general  form  of  the 
solution  exactly  satisfies  the  differential  equation,  the  bottom  boundary  condition,  and,  is 
properly  periodic  in  the  x-direction.  The  solution  does  not  provide  exact  fits  to  either  the 
kinematic  or  dynamic  free  surface  boundary  conditions.  Suppose  that  the  (n-l)fh  order 
solution  is  known  and  that  the  n^h  order  theory  is  to  be  developed.  The  nih  coefficients  arc 
determined  such  that  they  minimize  the  errors  in  the  two  free  surface  boundary  conditions 
at  the  (n-l)lh  order.  A significant  problem  is  that  the  configuration  of  the  n^h  order  water 
surface  is  not  known,  a priori;  it  is  therefore  necessary  to  best  satisfy  the  boundary 
conditions  on  an  approximate  expansion  of  the  nih  order  water  surface.  The  apparent  effect 
of  minimizing  the  errors  present  on  the  approximate  nl^  order  water  surface  is  that  the 
resulting  theory  of  a given  order,  if  convergent,  may  not  provide  the  best  fit  possible  for  the 
number  of  terms  (order)  included. 

As  a comparison  with  the  preceding  discussion  of  the  Stokes’  theory’,  consider  the 
corresponding  features  of  a Stream -function  theory  solution.  The  general  form  of  the 
solution  exactly  satisfies  all  of  the  boundary  value  problem  requirements  except  the 
DFSBC. 

At  this  stage,  one  inherent  advantage  of  the  Stream-function  theory  is  evident— all  of  the 
“free”  parameters  can  be  chosen  to  provide  a best  fit  to  the  DFSBC.  A second  and 
important  inherent  advantage  is  that  for  a given  nih  order  wave  theory,  all  of  the 
coefficients  are  chosen  such  that  they  best  satisfy  the  boundary  condition  on  the  nib  order 
water  surface.  The  distinction  is  that  because  a numerical  iteration  approach  is  used,  the  nib 
order  wave  form  is  known  (through  iteration)  at  that  order  of  solution.  Other  advantages  of 
the  Stream-function  theory  are  that  a solution  can  readily  be  obtained  to  any  reasonable 
order,  and  that  a measure  of  the  fit  to  the  one  remaining  boundary  condition  is  more  or  less 
automatically  obtained  in  the  course  of  the  solution.  Also,  the  form  of  the  terms  in  the 
solution  is  inherently  better  for  representing  nonlinear  waves,  due  to  the  t\  term  appearing  in 
the  argument  of  the  hyperbolic  sine  term  [cf.  Equation  (15)]. 

The  disadvantage  of  the  Stream-function  theory  is  that,  unless  tabulated  parameters  are 
available,  it  does  require  the  use  of  a digital  computer  with  a reasonably  large  memory.  The 
complexity  of  other  nonlinear  theories,  however,  generally  also  requires  the  use  of  a 
high-speed  computer. 

It  is  noted  that  a similar  but  different  Stream-function  theory  has  been  developed  and 
reported  by  Von  .Schwind  and  Reid  (1972)  subsequent  to  the  analytical  validity  study 
reported  here,  and  employs  a definition  of  the  DFSBC  error  which  differs  from  that  in  the 
present  study.  The  paper  by  Von  Schwind  and  Reid  presents  boundary  condition  errors  for 
three  wave  cases.  A comparison  between  their  errors  and  those  resulting  from  the 
Stream-function  theory  will  be  presented. 
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Analytical  Validity 

The  analytical  validity  of  a |)articiilar  wave  theory  has  In^en  previously  defined  as  the 
degre»‘  to  which  the  theory  satisfies  the  defining  equations,  i.e.,  Equations  (9)  through  (13). 
Again,  for  emphasis,  it  is  noted  that  a theory  providing  an  exact  fit  to  the  boundary 
conditions  would  have  a perfect  analytical  validity.  However,  due  to  a.ssumptions  of  ideal 
flow,  etc.,  in  the  formulation  of  the  problem,  a perfect  analytical  validity  does  not  ensure 
that  the  theory  would  provide  a good  re[)resenlation  of  laboratory  or  field  jihenomenon. 

The  rea.son  for  viewing  the  [)roblem  in  two  steps,  i.e.,  analytical  and  experimental 
vabdity,  is  that  the  results  of  the  analytical  validity  test  would  at  least  tend  to  indicate  the 
relative  applicability  of  the  available  wave  theories  for  particular  wave  conditions.  Also,  the 
results  would  provide  guidance  about  whether  the  most  fruitful  ap[)roach  would  be  directed 
toward  a more  representative  formulation  of  water-wave  theories  or  toward  the 
improvement  of  the  solutions  of  existing  formulations. 

[tefinition  of  Boundary  Condition  Errors 

Most  wave  theories  exactly  sati.sfy  the  governing  differential  equation  and  bottom 
boundary  condition,  although  some  of  the  solutions  only  approximately  satisfy  the 
(bfferential  equation.  Table  A lists  a number  of  the  theories  available  for  design  use  and  also 
indicates  the  conditions  of  the  formulation  which  are  satisfied  exactly  by  each  of  the 
theoiies.  Insjiection  of  Table  A shows  that  the  two  nonlinear  (free  surface)  boundary 
conditions  provide  the  Ix^st  basis  for  assessing  the  analytical  validity,  because  no  theory 
exactly  satisfies  both  of  these  conditions. 

Errors  based  on  the  dynamic  and  kinematic  free  surface  boundary  conditions,  are  defined 
as  functions  of  phase  angle  (0)  as  follows: 


El  (6) 


9n  _ w 
3x  u - C 


(16) 


£2(6)  = n + [ (u  - C)  ="  + w^']  - - Q (17) 

where  Q represents  the  mean  value  of  the  quantity  Q (Bernoulli  constant)  defined  in 
(Equation  (12).  Overall  errors  are  defined  as  the  root  mean  squares  of  the  distributed  errors. 


where  j represents  sampling  at  various  (evenly  spaced)  phase  angles. 


(18) 


(19) 
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TABLE  A 


Water  Wave  Theories  Included  in  Evaluation  Presented  by  Dean  (1968a) 


Theory 

Exactly  Satisfies 

DE 

BBC 

KFSBC 

DFSBC 

Linear  Wave  Theory— Airy 
(Ippen,  1966) 

X 

X 

— 

— 

Third  Order  Stokes 
(Skjelbreia  and  Hendrickson,  1961, 
as  summarized  by  Le  Mehaute  and  Webb,  1964) 

1 

X 

— 

Fiftli  Order  Stokes 
(Skjelbreia  and  Hendrickson,  1961) 

■ 

X 

— 

First  Order  Cnoidal 
(Laitone,  1960) 

■ 

X 

— 

— 

Second  Order  Cnoidal 
(Laitone,  1960) 

■ 

X 

— 

— 

First  Order  Solitary 

(Boussinesq,  as  summarized  by  Munk,  1949) 

X 

X 

— 

— 

Second  Order  Solitary 
(McGowan,  as  summarized  by  Munk,  1949) 

■ 

X 

X 

— 

Stream-Function  Numerical  Wave  Theory— Fifth  Order 
(Dean,  1968a) 

B 

X 

X 

— 

Results  of  Analytical  Validity  Comparison 

Most  of  the  results  of  the  study  of  analytical  validity  carried  out  under  this  project  have 
been  published  elsewhere  (Dean,  lQ6Ra),  and  therefore  will  he  reviewed  only  briefly  here. 

The  study  included  40  wave  cases  as  shown  in  Figure  2.  For  each  of  these  cases,  the 
overall  errors,  Ej  and  were  calculated  for  the  wave  theories  shown  in  Table  A.  The 
overall  dynamic  free  surface  boundary  condition  errors  were  made  dimensionless  by  dividing 
by  the  wave  height,  H,  i.e.. 


E2'  = Ei/H  (20) 

The  overall  kinematic  free  surface  boundary  condition  error  is  dimensionless  as  defined  in 
Equation  (18). 

Plots  of  the  dimensionless  kinematic  and  dynamic  free  surface  boundary  condition  errors 
arc  presented  in  Figures  3, 4,  5,  and  6 for  H/Hg  = 0.25  and  1.0  (Hg  = breaking  wave 
height).  The  KFSBC  error  is  identically  zero  for  the  Stream-function  and  McGowan  theories. 

As  stated  previously,  it  is  difficult  to  select  a single  index  that  would  clearly  be 
representative  of  the  overall  validity  of  all  wave  theories.  However,  an  index  was  chosen  that 
provided  a severe  test  for  the  Stream-function  theory,  and  yet  this  theory  provided  the  best 
general  analytical  validity. 

The  following  evaluation  plan  was  adopted,  the  results  of  which  would  be  somewhat 
biased  against  the  Stream-function  theory.  Most  of  the  wave  theories  do  not  satisfy  exactly 
either  the  DFSBC  or  KFSBC;  however,  the  Stream-function  theory  does  satisfy  exactly  the 
KFSBC.  It  therefore  seems  reasonable  that  if  the  Stream-function  theory  can  be  shown  to 
compare  favorably  against  other  theories  on  the  basis  of  only  the  DFSBC,  then  it  should 
provide  an  even  better  analytical  validty  than  the  comparison  shows. 

In  the  analytical  validity  investigation,  the  eight  wave  theories  in  Table  A were  examined. 
Because  the  fifth  order  was  the  highest  of  the  Stokes  theories  available,  it  was  arbitrarily 
decided  to  include  the  Stream-function  theory  only  to  the  fifth  order. 

The  evaluation  was  then  based  on  comparisons  presented  in  Figures  3, 4,  5,  and  6 and 
also  on  tbe  corresponding  figures  for  H/Hg  = 0.50  and  0.75,  which  are  not  presented  here. 
The  results  of  this  study  are  shown  in  Figures  7 and  8. 

Figure  7 presents  the  results  for  all  theories  excluding  the  Stream-function  theory.  It  is 
seen  that  the  Stokes  V theory  provide.8  the  best  fit  for  deep  water,  the  Airy  theory  provides 
the  best  fit  in  a part  of  the  transitional  and  shallow -water  ranges,  and  the  first-order  Cnoidal 
wave  theory  generally  provides  the  best  fit  in  the  shallow-water  range. 

Figure  8 presents  the  same  type  of  information.  Only  the  fifth-order  Stream -function 
theory  is  included  and  provides  the  best  fit  over  a wide  range  including  all  of  the  transitional 
and  deepwater  wave  regions  and  also  a significant  [)art  of  the  shallow-water  range  included 
in  the  comparison.  The  Airy  wave  theory  provides  the  best  fit  for  a small  part  of  the 
shallow-water,  near-breaking  waves  and  the  first-order  Cnoidal  wave  theory  provides  the  best 
fit  for  the  remainder  of  the  shallow-water  region. 
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h/Lo 


In  evaluating  tlie  resnltis  obtained  in  the  shallow-water  region,  it  is  noted  that  one 
eighth-order  Stream-function  theory  was  ealeulal«'<l  for  breaking  wave  conditions  and 
h/T^  = 0.1  foot/seeond^  as  shown  in  Figure  6.  This  figure  shows  that  the  use  of  higher  order 
Stream-function  theories  would  extend  the  range  of  best  validity  of  this  theory  to  shallower 
conditions  (Figure  8). 

Comparison  with  Stream-function  Theory  Developed  by  f'on  Schwind  and  Reid 
As  noted  earlier.  Von  Schwind  and  Reid  (1972)  have  developed  a Stream-function 
theory  with  basic  similarities  to  the  theory  employed  in  this  study.  The  principal  difference 
between  the  two  theories  is  that  Von  Schwind  and  Reid  transform  their  problem  to  and 
carry  out  their  solution  in  the  complex  plane.  It  is  noted  that  their  solution  in  terms  of 
wavelength  and  coefficients  is  also  obtained  by  iteration.  The  DFSIW.  error 
definition  e^(d),  used  by  V<m  Sehwind  and  Reid  was  originally  defined  by  (.happelear 
(1961),  and  is  somewhat  different  from  that  employed  here  (Equation  17)  and  is 


62  (0) 


£2  (0) 
Q + h 


(17a) 


It  is  noted  by  comparison  of  Equati<m8  ( 1 7)  and  (17a),  that  the  actual  distribution  of 
DFSBC  errors  would  appear  as  numerically  smaller  ha.sed  on  Equation  (17a)  due  to  the 
water  depth  and  Bernoulli  constant  appearing  in  the  denominator. 

Von  Schwind  and  Reid  presented  distributed  BFSBC  errors  for  three  sets  of  wave 
conditions.  Errors  were  calculated  for  the  same  wave  conditions  using  the  present  theory  . 
Figures  9,  10,  and  11,  are  reproduced  from  Von  Schwind  and  Reid,  and  the  maximum 
errors  obtained  by  the  present  theory  [indicated  University  of  Florida  (UF)]  are  shown  for 
each  wave  case.  The  maximum  UF  errors  obtained  are  so  small  that  it  would  not  he 
worthwhile  to  show  them  graphically.  Note  that  all  errors  (e^)  shown  in 
Figures  9,  10,  and  II  are  based  on  Equation  (17a).  The  reason  that  the  errors  obtained  by 
the  present  theory  are  smaller  than  those  obtained  by  Von  Schwind  and  Reid  is  not  known. 
With  a numerical  solution,  it  is  possible  to  obtain  a low  error  (down  to  some  limit)  by 
increasing  the  order  of  the  theory  or  by  increasing  the  number  of  iterations  used  to  obtain 
the  solution.  For  tlie  three  cases  shown  in  Figures  9 through  11,  the  UF  waves  were 
seventh-order  and  each  solution  was  obtained  by  15  iterations;  the  corresponding  values  for 
the  Von  Schwind-Reid  waves  are  not  known. 

Conclusions  Resulting  from  the  Analytical  Validity  Study 

Tl>e  analytical  validity  evaluation  is  based  on  the  degree  to  which  the  various  theories 
satisfy  the  governing  equations  in  the  boundary  value  problem  formulation.  It  is  stressed 
again  that  tliere  is  no  guarantee  that  a theory  providing  a good  analytical  validity  will 
necessarily  represent  well  the  features  of  the  natural  wave  phenomenon.  The  reason  is  that 
there  are  assumptions  (negligible  viscosity  and  capillary  effects)  introduced  into  the 
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Figure  9.  Comparison  of  errors  in  dynamic  free  surface  boundary  condition  for 
three  numerical  wave  theories,  wave  no.  1 
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Fifure  10.  Compariaon  of  errors  in  dynamic  free  surface  boundary  condition  for 
three  numerical  wave  theories,  wave  no.  2 
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governing  equations  which  may  adversely  affect  the  degree  to  which  the  formulation 
represents  real  wave  motion.  The  purpose  of  the  analytical  validity  study,  rather,  was  to 
attempt  to  resolve  the  question  of  whether  the  theories  developed  for  the  same  formulation 
and  for  various  regions  of  relative  depth  do  indeed  provide  the  best  fit  in  these  regions.  Also 
this  study,  combined  with  some  additional  studies  reported  later  in  this  report,  does  aid  in 
determining  whether  the  most  critical  need  in  wave  theory  research  is  in  the  improvement  of 
the  formulation  or  in  the  development  of  improved  solutions  to  the  existing  formulation. 

The  results  of  the  analytical  validity  study  have  shown  that: 

1.  The  general  status  of  wave  theories  for  h/t^  > 0.2  foot/second* , for  instance,  is  much 

more  satisfactory  than  for  the  smaller  values  of  h/T^ . In  particular,  for  the  larger  relative 
depths,  there  is  reasonable  consistency  between  the  fits  to  the  dynamic  free  surface 
boundary  condition  and  the  maximum  drag  force  as  calculated  by  the  various  theories 

including  a seventh-order  Stream-function  theory.  In  shallow  water,  it  is  not  clear  that  the 

boundary  condition  fit  is  an  appropriate  measure  of  wave  theory  validity,  unless  the 

associated  errors  are  very  small.  In  particular,  the  Airy  wave  theory  provides  a relatively 

good  fit  to  the  boundary  conditions  in  shallow  water;  however  this  theory  does  not 
represent  many  of  the  observed  features  of  shallow-water  waves  including  the  strong 
skewness  of  the  wave  profile  about  the  mean  water  level. 

2.  The  Stokes  higher  order  wave  theories  converge  to  accurate  representations  of  wave 
motion  in  deep  water;  however,  in  transitional  and  shallow  water,  the  boundary  condition 
fits  are  relatively  poor.  Furthermore,  no  fifth-order  Stokes  theory  solution  could  be  found 
for  shallow-water  waves  or  the  smaller  values  of  the  transitional  zone.  The  limiting  value 
of  h/T*  for  which  a solution  exists,  depends  on  H/T*  and  was  in  the  range  of 
0.1  < h/T*  < 0.5  foot/second*  for  the  conditions  examined. 

3.  Finally,  it  is  observed  that  the  second-order  Cnoidal  theory  provided  a worse  fit  to  the 
boundary  conditions  than  the  Hrst -order  Cnoidal  theory  for  all  wave  conditions  examined. 
There  are  other  versions  of  Cnoidal  theories;  the  boundary  condition  fits  of  these  theories 
have  not  been  evaluated  in  this  study. 

4.  The  Stream-function  theory  described  in  this  report  provides  good  analytical  validity 
over  a wide  range  of  wave  conditions. 

The  reader  is  referred  to  Dean,  (1968a)  for  reinforcement  of  statements  presented. 

Experimental  yalidity 

As  previously  described,  experimental  validity  is  based  on  the  comparison  of  theoretical 
predictions  and  measured  wave  phenomena.  If  it  could  be  generally  shown  that  the  theory 
providing  the  best  analytical  validity  also  provides  the  best  experimental  validity,  then  it 
could  be  concluded  that  the  formulation  is  valid  and  that  the  errors  in  the  boundary 
conditions  are  also  good  indicators  of  experimental  validity.  If  the  differences  between  the 
theory  and  experiments  were  of  the  same  order  as  the  estimated  experimental  error,  and  if 
this  could  be  shown  to  be  the  situation  generaUy,  then  the  most  productive  direction  in 


22 


water-wave  research  on  this  problem  would  be  improved  measurements.  If  however,  the 
disagreement  ludween  theory  and  experiment  is  much  larger  than  can  be  attributed  to 
ex|)erimental  error,  and  especially  if  this  difference  were  of  engineering  significance,  then 
additional  efforts  on  the  formulation  and  solution  of  water  wave  theories  would  be 
indicated. 

The  availability  of  data  is  inadecpiate  to  carrv  out  a comprehensive  evaluation  of 
exjHTimental  validity  over  all  ranges  of  relative  depth  and  heights  of  engineering  importance, 
be  Mehaute,  Divorky,  and  Lin  (1968)  have  carried  out  a measurement  program  in  which 
distributions  over  depth  of  horizontal  water  particle  velocities  were  measured  under  the 
crest  phase  position  of  fairly  high  waves  in  the  .shallow  and  transitional  depth  ranges.  The 
residts  included  measured  horizontal  water  particle  velocity  distributions  for  eight  cases,  and 
also  a vertical  water  particle  velocity  distribution  for  one  case,  and  one  measured  wave 
profile.  I.e  Mehaute,  Divorky,  and  Lin  compared  a number  of  wave  theories  with  their  data; 
however  the  .Stream-function  theory  was  not  included.  The  ex}*erimental  validity  reported 
in  this  study  was  based  on  a comparison  of  the  Stream-function  theory  with  the  data 
described  earlier. 

It  should  be  emphasized  that  the  only  addition  to  the  paper  by  Le  Mehaute,  Divorky,  and 
Lin  (1968)  is  ( 1 ) comparison  of  the  Stream-function  theory  with  the  data  and 
(2)  calculations  wluch  repre.sent  the  overall  agreement  between  the  data  and  several  of  the 
theories.  In  the  Stream-function  horizontal  velocity  component  profiles  presented,  a 
uniform  mass  transj)ort  velocity  has  been  subtracted  out,  whereas  due  to  time  limitations, 
the  other  theoretical  velocity  distributions  were  .simply  plotted  from  Le  Mehaute,  Divorky, 
and  Lin.  It  is  not  clear  whether  or  not  the  mass  transport  term  should  be  subtracted  out. 
Although  the  experiments  were  conducted  in  a closed  tank,  the  data  were  taken  before 
waves  reflected  from  the  beach  had  propagated  back  to  the  tank  test  section,  and  the  zero 
net  flow  over  depth  had  probably  not  been  establi.shed  completely. 

In  all,  data  for  10  different  wave  conditions  are  available.  These  waves  are  in  the  shallow 
and  transitional  dejith  regions,  and  according  to  the  conventional  breaking  criteria,  the  wave 
heights  range  from  0.43  to  0.70  of  the  breaking  height.  The  wave  conditions  are  shown  as 
points  in  Figure  12  where  isolines  representing  various  ratios  of  wave  height  to  breaking 
wave  height  are  also  presented.  It  is  emphasized  that  the  breaking  wave  height  in  Figure  12 
is  the  conventional  breaking  height:  i.e.,  H/h  = 0.78  in  shallow  water  (McGowan  reviewed  by 
Munk,  1949);  H/L  = 0.142  in  deep  water  (Michell,  1893);  in  the  transitional  range,  the 
breaking  limit  was  first  establi.slied  by  Reid  and  Bretschneider  (1953)  by  interpolating  on 
the  basis  of  measured  data  and  is  presented  in  several  more  available  references,  e.g.  (Ippen, 
1966)  and  (Bretschneider,  1960).  A recent  pa[)er  by  Divoky,  Le  Mehaute,  and  Lin  (1970) 
reports  an  experimentally  determined  shallow-water  breaking  limit  of  approximately 
Hg/h  = 0.60  to  0.66  as  compared  to  the  conventional  value  of  0.78.  The  recent  experiments 
resulting  in  the  lower  value  were  obtained  with  a laterally  converging  wave  channel. 
Certainly  it  is  apparent  that  more  work  is  needed  to  better  resolve  wave  breaking  limits. 
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Tabic  B presents  the  comparison  results  included  in  the  experimental  validity  evaluation. 
The  eight  comparisons  of  horizontal  water  particle  velocity  are  presented  in  Figures  13 
through  20;  the  vertical  velocity  comparison  is  presented  in  Figure  21;  and  the  wave  profile 
is  presented  in  Figure  22. 

Figures  1 3 through  20  indicate  that  the  Stream-function  theory  is  in  reasonable 
ap'eement  with  the  data.  It  is  noteworthy  that  the  shallow-water  wave  theories  which 
should  provide  good  fits  to  the  data  are  so  poor.  Another  interesting  feature  of  the 
comparison  is  that  the  linear  (Airy)  wave  theory  agrees  better  with  the  data  tiian  would  be 
expected. 

Of  the  12  theories  included  in  the  comparison,  the  better  agreements  with  data  were 
provided  by  the  following  five  theories:  Airy,  Keulegan  and  Patterson  Cnoidal  wave  theory, 
Goda,  Long-Wave,  and  Stream-function.  These  five  theories  were  then  selected  for  further 
examination  of  their  agreement  with  the  data.  The  standard  deviations  between  each  of 
these  theories  and  the  data  were  calculated,  and  are  presented  in  Table  C where  it  is  seen 
that  tlie  Stream-function  theory  provided  the  best  fit  to  the  data,  followed,  in  order,  by  the 
Goda,  Keulegan  and  Patterson  Cnoidal,  Airy,  and  the  Long-Wave  Theories. 

The  Goda  “theory”  is  actually  a series  representation  in  which  the  analytical  forms  of  the 
terms  comprising  the  series  are  tlie  same  as  the  hyperbolic  and  trigonometric  functions  in 
the  Stokes  theories.  However,  the  coefficients  modifying  these  terms  were  determined 
empirically  by  wave  tank  experiments. 

Additional  calculations  not  presented  here  showed  that,  assuming  the  data  were  valid,  the 
Stream-function  wave  theory  would  on  the  average  overprcdict  the  maximum  total  drag 
force  on  a vertical  cylinder  by  21  percent. 

Data  representing  the  vertical  velocity  distribution  with  depth  are  available  for  only  one 
set  of  wave  conditions,  (see  Figure  21).  The  McGowan  theory  provides  the  best  fit  to  the 
data;  the  next  best  fit  is  associated  with  the  .Stream-function  theory.  Differences  between 
the  McGowan  and  Stream-function  theories,  however,  are  quite  small  and  it  is  probably  not 
justified  to  draw  conclusions  from  only  one  set  of  data.  Interpreted  in  terms  of  vertical  drag 
forces  on  a horizontal  cylinder,  the  Stream-function  would  underpredict  the  forces  by  30 
percent. 

Tlie  one  set  of  wave  profile  data  are  compared  with  the  various  theories  in  Figure  22. 
Although  no  detailed  comparisons  were  made,  it  appears  that  the  Stream -function  theory  is 
in  as  good  or  better  agreement  than  any  of  the  other  theories  shown. 

Conclusions  Resulting  from  the  Experimental  Validity  Study 

Gomparisons  of  Stream-function  theory  predictions  with  measurements  of  velocity 
components  and  one  wave  form  representing  transitional  and  shallow-water  waves  indicate 
reasonably  good  agreement.  Interpreted  on  the  ba.sis  of  maximum  horizontal  drag  force 
components,  the  Stream-function  theory  would  over  predict  by  an  average  of  21  percent. 
Recognizing  that  the  experimental  accuracy  is  approximately  5 percent  these  results  are 
considered  reasonable  for  engineering  applications.  The  predicted  maximum  vertical  drag 
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Horizontal  water  particle  velocity  under  the  crest,  Case  2 


lonzontal  water  particle  velocity  under  the  cresi 


Figure  16.  Horizontal  water  particle  velocity  under  the  crest.  Case  4 
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Figure  18.  Horizontal  water  particle  velocity  under  the  crest.  Case  6 


Figure  20.  Horizontal  water  particle  velocity  under  the  crest.  Case  8 
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Figure  21.  Vertical  water  particle  velocity,  Case  9 
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DISTANCE  FROM  CREST,  x (ft) 


TABLE  C 


Standard  Deviation  of  Differences  Between  Horizontal 
Velocities:  Measured  vs.  Predicted 


Case  No. 

Standard 

Deviation , 

a (ft/sec) 

Theory 



Airy 

Long  Wave 

God  a 

K & P Cnoidal 

1 

0.229 

0.232 

0.328 

0.413 

0.396 

2 

0.139 

0.234 

0.297 

0.146 

0.211 

3 

0.096 

0.470 

0.468 

0.206 

0.155 

4 

0.126 

0.442 

0.453 

0.134 

0.136 

5 

0.245 

0.225 

0.291 

0.357 

0.487 

6 

0.216 

0.181 

0.244 

0.095 

0.469 

7 

0.123 

0.493 

0.513 

0.316 

0.188 

8 

0.183 

0.418 

0.434 

0.215 

0.272 

Average 

0.170 

0.337 

0.379 

0.235 

0.289 

u„  = measured  velocity  component 
M 

u^  ® theoretical  velocity  component 
J = number  of  levels  considered  for  each  case  (14  to  15) 


i 
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forces  on  a horizontal  cylinder  would  he  loo  small  li\  .'50  ixtcciiI:  Iiowcvct,  this  statement  is 
based  on  a comparison  w^ith  only  one  set  of  data.  Good  (jualitative  agreement  was  found 
between  measured  and  predicted  wave  profiles  • 

Finally,  based  on  the  results  of  both  the  analytical  and  experimental  validity  studies,  it  is 
concluded  that  the  Stream-function  theory  is  best  suited  for  engineering  design  purposes.  It 
was  decided  to  tabulate  variables  that  would  be  of  use  in  engineering  design  as  calculated 
from  the  Stream-function  theory.  The  next  section  describes  the  variables  included  in  the 
tables. 


IV.  DESCRIPTION  OF  TABLES 

Introduction 

An  attempt  has  been  made  to  include  in  the  tables  those  variables  of  greatest  present 
engineering  interest  and  application.  In  addition,  other  variables  were  included  which  would 
be  relevant  to  checking  the  relative  analytical  validity  of  other  theories  or  variables  which 
were  of  scientific  interest  and  could  conceivably  be  required  for  engineering  in  the  future. 
V ariables  have  been  included  which  describe  the  detailed  kinematics  of  the  waves  and  also 
wliich  represent,  e.g.,  the  integrated  effect  of  the  flow  on  a structural  member. 

It  is  not  possible  to  assemble  in  concise  tabular  form  all  variables  that  could  be  of 
engineering  use.  It  is  feasible  to  tabulate  the  dimensionless  drag  force  for  all  vertical  piling 
extending  from  the  bottom  up  to  a certain  level.  It  would  not  be  feasible,  however,  to 
concisely  tabulate  the  total  drag  force  on  members  with  all  possible  inclinations  relative  to  a 
vertical. 

Forty  sets  of  dimensionless  wave  conditions  were  selected  for  tabulation.  Each  case  is 
characterized  by  values  of  h/L^  and  H/L^,.  The  parameter  h/L^,  ranged  from  0.002  to  2.0 
and  covered  the  relative  depth  range  from  shallow  to  deep  water.  The 
parameter  H/L^  included  wave  steepnesses  ratios:  0.25,  0.5,  0.75,  and  1.0  of  the  breaking 
wave  steepness  for  each  of  the  10  h/L^  values  tabulated.  Figure  23  shows  the 
dimensionless  wave  conditions  selected  for  tabulation  and  also  indicates  the  referencing 
notation  for  the  cases. 

All  tabulated  variables  are  presented  in  dimensionless  form.  The  description  of  these 
variables  is  presented  in  the  following  paragraphs  and  in  Tables  D,  E,  and  F,  where  generally 
the  following  are  included:  the  equation  for  the  variable,  the  dimensionless  form  of  the 
variable,  an  equation  number  for  reference  purposes,  and  tlie  table  number  in  the  wave 
tables.  To  reduce  confusion,  it  should  be  noted  that  the  tables  presented  in  this  report  are 
denoted  by  letters;  the  wave  tables  are  identified  by  Roman  numerals. 

Variables  Presented  in  Tabular  Form 

Three  classes  of  variables  are  tabulated:  (1)  Internal  field  variables,  depending  on  6 and  S, 
(2)  Variables  depending  on  6 only,  and  (3)  Overall  variables  which  have  a single  value  for 
the  entire  wave,  for  example  the  wavelength. 
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Figure  23.  Wave  characteristics  selected  for  tabulation 
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Inlernat  Field  I'arinblen  l)ei>eitdinp  on  0 and  S 

1'h<‘  inirriial  field  variables  are  talinlaled  at  eijually  spac'ed  dimensionless  distances  above 
the  bottom,  i.e.,  at  .''/b  values  of  0,0. 1,0.2  . . . uj)  to  and  including  the  free  surface,  and 
at  0 values  of  0°.  10°,  20°,  30°,  50°,  75°,  100°,  130°,  1B0°.  Figure  24  shows  a .sample 
pre.sentatioii  of  the  dimensionless  horizontal  velocity  component  field. 

.‘\  description  of  the  entries  in  Figure  24  will  s«’rve  to  familiarize  the  reader  with  most  of 
the  features  of  the  tables.  The  first  row  lists  the  phase  angles  (theta)  in  degrees.  The  second 
row  lists  the  dimen.siunle.ss  wave  profile  (tj/H)  at  the  corresponding  pha.se  angles.  The 
|)ereent  values  listed  heneath  the  tj/H  values  are  the  differences  between  the 
.‘stream-function  and  Airy  Theories,  defined  as: 


Percent 


Stream-function— Airv 

X 100  percent 

Stream-function 


The  main  body  (remaining  portion)  of  the  table  lists  the  dimensionless  horizontal  water 
particle  velocities.  The  row  labeled  “.Surface”  represents  the  dimensionless  velocities 
evaluated  at  the  free  surface;  the  jierccritagc  differences  for  velocities  are  calculated  as 
defined  above  for  the  profile.  The  remaining  part  of  the  tabh’  r<-pres<*nts  the  dimensionle.ss 
velocities  and  perirntage  differences  evaluated  on  a grid  of  (0,  S/h),  The  lack  of  entries  for 
the  higher  .S/h  and  higher  tlieta  values  (right  side  of  page)  results  from  the  wave  profile  in 
the  trough  region  lieing  lower  than  in  the  cre.st  region  (left  side  of  page).  Two  additional 
comments  pertaining  to  the  percentage  values  will  complete  the  description  of  the  sample 
table.  ,A  jieri’entage  difference  value  of  exaetly  100  jK-reent  implies  that  the 
Stream-function  profile  occured  at  a (0,S/h)  value,  however,  the  Airy  profile  was  lower 
than  the  particular  S/h  at  the  phase  angle,  0,  i.e.,  this  grid  point  was  not  “covered” 
by  the  Airy  pr«)file.  For  example,  this  is  the  ca,se  at  0 = 0°,  S/h  = 1 .5  and  1.6  and  0 = 180°, 
.S/h  = O.H  and  0.9.  F'inally,  the  asterisks  indicate  that  the  percentage  differences  were  not 
calculated  beeau.s<“  the  Stream-function  value  was  le.ss  than  5 percent  of  the  maximum 
.Stream-function  value.  This  avoided  the  tabulation  of  very  large  j)erecntages  which  would 
have  been  the  result  of  division  by  a small  number. 

A brief  description  of  each  of  the  tabulated  internal  field  variables  is  presimted  below. 

Horizontal  Water  Particle  Velocity  Comitonent,  u{6 ,S) 

The  horizontal  water  particle  vek>eity  component,  u(0,  S),  is  defined  by  Equation  (21). 
(The  equations  for  the  tabulated  functions  are  presented  in  Tables  1),  E,  and  F.)  The  values 
u'(0,  S)  tabulated,  are  presented  (Table  I)  in  the  following  dimensionless  form: 


dimensionless  horizontal  velocity  component  field 


Vertical  Water  Particle  Velocity  Comfjonent,  w(0,  S) 

The  vertical  water  particle  velocity  component,  w(0,  S),  is  defined  by  Equation  (22). 
The  dimensionless  values  tabulated  (Table  II),  w'(0,  S),  are  defined  by: 


w'  (6,  S) 


w(9,  S) 
(H/T) 


Horizontal  Water  Particle  Acceleration,  Du/Dt 

The  horizontal  water  particle  acceleration,  Du/Dt'  is  defined  in  terms  of  the  velocity 
components  as  presented  in  Equation  (23).  Note  that  the  tabulated  values  represent  the 
total  (or  material,  substantial,  etc.)  acceleration  consisting  of  the  sum  of  the  local  and 
advective  contributions.  The  dimensionless  values  tabulated  (Table  III),  Du'/Dt',  are 
defined  by: 


Du ' _ 1 Du 

Dt  " (H/T^)  Dt 

Vertical  Water  Particle  Acceleration,  Dw/Dt 

The  vertical  water  particle  acceleration,  defined  in  Equation  (24),  is  tabulated 
(Table  IV)  in  the  following  dimensionless  form: 


Dw^  _ 1 Dw 

Dt  (H/T^)  Dt 


Drag  Force  Component,  Fq{6,  S) 

The  drag  force  component  up  to  a certain  elevation,  S,  is  defined  by  Equation  (25) 
and  tabulated  (Table  V)  in  dimensionless  form  as: 


^D 


C^pD(H/T) 


Fd 


Inertia  Force  Component,  Fj(B,  S) 

TTie  inertia  force  component  up  to  a certain  elevation,  S,  is  defined  by  Equation  (26) 
and  tabulated  (Table  VI)  in  dimensionless  form  as: 


C^pTiD*  (H/T*)h 
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Drofi  Moment  ('.oniftonent.  MpiO,  S) 

TIu'  drag  nmment  coinpotifnl  about  llir  bottom  du«;  to  wave  pressures  acting  on  a 
vertical  member  extending  up  to  an  elevation,  S,  is  j)resented  as  F^quation  (27)  and 
presented  (Table  VTI)  in  dimensionless  form  as: 


CdPD(H/T)  "h" 


Mr 


Inertia  Moment  Component.  MjlO,  S) 

The  inertia  moment  component  about  the  bottom  due  to  wave  pressures  acting  on  a 
vertical  member  extending  up  to  an  elevation  S,  is  defined  in  Equation,  (28)  and  |)re8ented 
(Table  VTIl)  in  dimensionle.ss  form  as; 


m;  = 


Cj^PTtD"  (H/T^h^ 


Ml 


Ilynamie  f*ressure  Component, 

The  dynamic  pressure  component,  defined  by  Equation,  (29)  is  tabulated  (Table  IX)  in 
dimensionless  form  as: 


Pd 


f 


Pd 


Thts  completes  the  description  of  the  field  variables  (depending  on  6 and  S)  that  are 
included  in  the  tables. 

Variables  Depending  on  6 Only 
Hater  Surface  Displacement,  rj{0) 

The  water  surface  displacement  is  defined  in  Equation  (30),  and  tabulated 
(Tables  I through  IX)  in  dimensionless  form  as: 


Total  Drag  Force  Component,  Fp{9) 

The  total  drag  force  component  is  defined  by  E(juation  (25)  with  the  upfter  limit  taken 
to  be  h + r?(0),  and  is  tabulated  (Table  V,  labeled  “SURFACE”)  in  dimensionless  form  as: 


Fd'  = 


CdPD(H/T) ^hj 


Fd 
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Total  Inertia  Force  Component,  Ff{0) 

The  total  inertia  force  component  is  defined  by  Equation  (26)  with  the  upper  limit 
taken  to  be  h + r)(6),  and  is  tabulated  (Table  VI,  labeled  "SURFACE”)  in  dimensionless 
form  as: 


C^pTrO^  (H/T2)hJ 


Fl 


Total  Drag  Moment  Component,  Mfy(d) 

The  total  drag  moment  component  is  defined  by  Equation  (27)  with  the  upper  limit 
taken  to  be  h + tj(6)  and  is  tabulated  (Table  VII,  labeled  "SURFACE”)  in  dimensionless 
form  as: 


”d'  = 


C^pD (H/T) 


Mr 


Total  Inertia  Moment  Component,  Mf(6) 

The  total  inertia  moment  component  is  defined  by  Equation  (28)  with  an  upper  limit 
of  b + T](d)  and  is  tabulated  (Table  VIII,  labeled  "SURFACE”)  in  dimensionless  for  as: 


Mi'  = 


4 

CMPtrO^"  {H/T2)h" 


Mi 


Kinematic  Free  Surface  Boundary  Condition  Error,  e,  (0) 

The  kinematic  free  surface  boundary  condition  error  is  defined  by  Equation  (35).  This 
variable,  as  defined,  is  in  dimensionless  form  and  is  tabulated  in  Table  \: 

Item  1 , Linear  Wave  Theory- 
Item  2,  Stream-Function  Theory 

Dynamic  Free  Surface  Boundary  Condition  Error,  6^(0) 

The  dynamic  free  surface  boundary  condition  error  is  defined  by  Equation  (36)  and  is 
tabulated  (Table  X)  in  the  following  dimensionless  form: 


£2  ' = 


G2 


with 

Item  3,  Linear  Wave  Theory 
Item  4,  Stream-Function  Theory 

Ihis  completes  the  presentation  of  variables  depending  on  0 only. 
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Overall  Variables  (do  not  depend  on  0 or  S) 

Wavelength,  L 

For  the  Stream-function  theory,  there  is  no  definable  expression  for  the  wavelength. 
Rather  the  wavelength  is  determined  as  a part  of  the  numerical  solution  as  described  in 
Appendix  I.  The  dimensionless  wavelength  is  presented  (Table  XI,  Item  1)  in  the  following 
dimensionless  form: 


L'  = 


L 


Average  Potential  Energy,  PE 

The  average  potential  energy  is  defined  by  Equation  (38)  and  is  tabulated  (Table  XI, 
Item  2)  in  dimensionless  form  as: 


PE' 


PE 


Note  that  the  dimensionless  form  is  defined  to  be  0.5  for  the  linear  (Airy)  wave  theory. 
Average  Kinetic  Energy,  KE 

The  average  kinetic  energy  is  defined  by  Equation  (39),  and  b also  tabulated  (Table  XI, 
Item  3)  in  dimensionless  form  as: 


KE' 


KE 


As  for  the  dimensionless  potential  energy,  the  dimensionless  value  for  the  linear  (Airy) 
wave  theory  is  0.5. 

Average  Total  Energy,  TE 

The  average  total  energy  is  simply  the  sum  of  the  potential  and  kinetic  energy 
contributions  (Equation  (40)],  and  is  tabulated  in  dimensionless  form  (Table  XI,  Item  4) 
such  that  the  difference  from  unity  is  an  indication  of  the  deviation  from  the  linear  wave 
theory. 


TE' 


Average  Total  Energy  Flux,  Fj£ 

Tlte  average  total  energy  flux  is  defined  by  Equation  (41),  and  is  tabulated  (Table  XI, 
Item  5)  in  dimensionless  form  as: 


8 

YH='  l/t 


TE 


50 


Group  VelocUy,  Cq 

Tlie  group  velocity  is  defined  as  the  ratio  of  total  energy  flux  to  total  energy 
[Equation  (42)]  and  is  presented  (Table  XI,  Item  6)  in  dimensionless  form  as: 

“ (fTr) 

The  dimensionless  group  velocity  is  defined  such  that  for  linear  wave  theory  the  shallow  and 
deepwater  values  are  1 .0  and  0.5,  respectively. 

Average  Momentum,  M 

The  total  average  momentum  is  defined  by  Equation  (43)  and  is  presented  (Table  XI, 
Item  7)  in  dimensionless  form  as: 


V.  J 


Tlie  dimensionless  momentum  is  defined  such  that  for  linear  wave  theory  the  result  is  unity. 
Note  that  mass  transport  velocity,  U = [M/ph)  is  proportional  to  the  average  momentum. 
Average  Momentum  Flux  in  fFave  Direction, 

The  total  average  momentum  flux  in  the  wave  direction  is  defined  by  Equation  (44) 
and  is  tabulated  (Table  XI,  Item  8)  in  the  following  dimensionless  form: 


The  above  definition  reduces  to  1.5  and  0.5  for  linear  wave  theory  for  shallow  and 
deepwater  waves,  respectively. 

Average  Momentum  Flux  Transverse  to  Wove  Direction,  F^ 

The  total  average  momentum  flux  in  a direction  perpendicular  to  the  wave  advance 
direction  is  needed  to  define  the  radiation  stress  tensor,  (discussed  by  Bowen  1969),  is 
defined  by  Equation  (45),  and  is  tabulated  (Table  XI,  Item  9)  in  the  following  dimensionless 
form: 


F' 

m 


y 


For  linear  wave  theory,  the  above  deTinition  reduces  to  0.5  and  0.0  for  shallow  and 
deepwater  waves,  respectively. 
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Kini'ruulii  Frri'  Siirfuc)'  lioundnry  ('.ondition  Frrors,  Cj 

'I'lic  kiiii-iiiiilic  frrc  siirl'acf  hoiitulary  condition  iTfor  is  dcfiticd  in  dimensionless  form 
li\  Ktjnalion  (do)  und  the  root  incan-s(|uarc  (KMS)  and  maximum  values  are  tabulated 
4 ( Table  XI.  Items  It)  and  12)  as  defined  by  Equation  (46). 

Itynamir  Frrv  Siirfarv  Htnindary  Ciniditinn  Frr»rs. 

'Tin*  dMiamie  free  surfaee  lioundarv  condition  error  is  defined  by  E(|uation  (.46)  and  is 
re|ires<-nled  in  tlie  followin;;  ditnensiordess  fortri: 


^2 

= ir 

The  HM.'s  and  maximum  values  are  tabulated  (Tabl(‘  XI,  Items  II  and  1,4)  as  defined  by 
K(|uatioti  (47). 

Kinonuttu'  Frrr  Surface  Itrcakirifi  Parameter, 

The  kinematic  free  surface  breakino;  parameter  is  tabulated  (Table  XJ,  Item  14)  as 
defined  by  Ecpiation  (4H)  (dimensiotde.ss  form). 

Itynamir  Free  Surface  Hreaking  Parameter,  (3^ 

The  dynamic  free  surfaee  breaking  parameter  is  tabidated  (Table  XI,  Item  1.4)  as 
defined  by  Fapjation  (40)  in  dimensionless  form. 

Variables  Presented  in  Gra/diical  Form— Combined  Effect  of  Shoaling  and  Refraction 

In  addition  to  developing  the  tabulated  values  previously  described,  the  study  included 
tile  development  of  the  combined  effect  of  shoaling  and  refraction  for  nonlinear  waves 
advancing  toward  shore  with  a deepwater  direction,  a„,  over  bathymetry  characterized  by 
straight  and  parallel  contours, 

E’or  linear  wave  theory,  it  is  po.ssible  to  .separate  the  .shoaling  and  refraction  effects, 
lieeause  neither  wave  celerity,  ('.  (governing  refraction),  nor  group  velocity,  Cq  (governing 
energy  flux),  is  dc|»endent  on  wave  height.  E’or  nonlinear  waves,  both  celerity  and  group 
velocity  at  a certain  location  depend  on  wave  height  as  well  as  wave  period  and  water  depth. 
The  shoaling-refraetion  effects  for  nonlinear  waves  are  therefore  not  .separable,  and  the 
combined  effect  depi-nils  on  the  deepwater  wave  steepne.ss,  as  well  as  the  local 

relative  dejith. 

Ib’causi'  the  .shoaling-refraetion  re.sults  are  not  readily  jiresented  in  tabular  form,  graphs 
are  pres«uited  as  Figures  2.4,  through  29  for  deepwater  wave  directions,  a„  of  0°,  10°,  20° 
40°,  and  60°.  A brief  des<Tiption  of  the  use  of  thes*'  graphs  follows.  A wave  with  a 
deepwater  dire<'ti(»n  a„,  will  propagate  toward  .shore  such  that  the  local  H/L^  will  fall 
along  a curve  characterized  by  the  deepwater  value  H^/L^.  At  any  particular  relative 
depth,  h/b^j,  the  local  wave  stec[iness  ll/b„  and  direction  a arc  read  from  the  ordinate 
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Combined  shoaling-refraction  for  a dcc{)water  wave  direction. 
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Figure  29.  Combined  shoaling-refraction  for  a deepwater  wave  direction, 


and  interpolated  from  the  appropriate  isolines,  respectively.  The  region  to  the  lower  right  of 
the  Une  of  dots  indicates  the  region  where  use  of  the  linear  theory  agrees  with  the  nonlinear 
results  presented  within  1 percent  in  H/L^  and  1°  in  wave  direction,  a. 

V.  EXAMPLES  ILLUSTRATING  USE  OF  WAVE  TABLES 

Introduction 

Tlie  preceding  chapter  has  described  the  formats  and  the  various  dimensionless 
parameters  included  in  the  wave  tables.  To  aid  in  the  application  of  the  tables,  examples  will 
be  presented  illustrating  their  use.  The  first  example  is  a [irohlem  of  a near -breaking  wave 
interacting  with  an  offshore  structure  supported  by  cylindrical  [tiling.  This  example  will  use 
those  tables  which  contain  the  wave  profile  and  the  wave  forces  and  moments.  Additional 
examples  will  then  be  presented  which  will  illustrate  the  use  of  most  of  the  remaining  wave 
tables.  Where  possible,  examples  were  selected  to  parallel  |)roblems  which  may  occur  in 
offshore  design. 

It  is  worthy  of  note  that  the  tables  have  a much  wider  ap|)licability  than  can  be  illustrated 
by  the  limited  number  of  examples  presented  here.  A thorough  familiarity  with  the 
information  summarized  in  the  tables  should  aid  in  an  understanding  of  them  and  their  use 
ill  many  problems  involving  water-wave  phenomena.  Tlie  examples  will  be  presented  in  U.S. 
Customary  units;  however  the  tables  are  in  dimensionless  form,  and  any  system  could  be 
used  readily. 

Example  1— Deck  Elevation  and  Wave  Forces  and  Moments  on  an  Offshore  Platform 

(Consider  the  design  problem  of  determining  the  deck  elevation  and  horizontal  wave 
forces  and  moments  upon  individual  members  of  the  offshore  platform  illustrated  in 
Figure  .30.  Suppose  that  the  design  depth  (mean  low  water  + maximum  tide  + storm 
surge),  h,  is  41  feet,  and  the  main  structural  members  of  the  jilatform  and  outriggers  are 
piUngs  6 feet  in  diameter,  wit*"  filing  fenders  3 feet  in  diameter.  The  fenders  extend  from 
4.1  feet  above  the  -design  Stillwater  level  to  a depth  of  8.2  feet.  The  outriggers  are  20.5  feet 
high.  Suppose  that  analysis  indicates  that  the  design  wave  will  have  a (breaking) 
height,  H.  of  31.78  feet  and  a jieriod,  T,  of  20  seconds.  The  drag  and  inertia 
coefficients,  Cjp  and  Cj^,  for  this  structure  are  assumed  to  be  1.05  and  1.5,  respectively. 

To  determine  which  set  of  tables  to  use,  calculate  h/L^  and  H/L^j,  where  =gT^/(2ir), 

Lo  (5.i2)  (20)' 

E?  ■ (S.li)  (20)^  - 


In  this  and  most  subsequent  examples  in  this  chapter,  the  tables  for  Case  4-D  will  be  used 
(see  Figure  23).  A sample  table  set  for  Case  4-D  is  included  as  Appendix  III. 


► 


1 

, /• 


Deck  Elevation 

To  fusun*  that  the  deck  is  above  tlie  design  crest  elevation;  thereby  avoiding  unnecessarily 
large  horizontal  and  vertical  forces  and  damage  to  the  platform  has*-,  the  height  of  the  lower 
elevation  of  the  deck  will  b«‘: 


h"  = h + 


n 

max 


+ h' 


where  h is  the  design  water  depth.  is  the  maximum  displacement  of  the  wave  above 

design  Stillwater  level,  and  h'  is  the  deck  freeboard  (say  10  feet  for  this  problem).  Vnax 
will  occur  at  zero  phase  angle  (0  = 0°)  and  from  any  of  the  first  nine  tables,  eta/height  = 
0.89  for  0 = 0°.  Therefore,  = 0.89  (H)  =28.3  feet  and  h"  = h + rj  + h' = 41  + 28.3 + 
10  = 79.3  feet.  The  platform  will  Im‘  constructed  so  the  lower  deck  elevation  will  be  79.3 
feet  above  the  bottom. 

In  determining  the  forces  and  moments,  it  is  assumed  that  the  piling  arc  sufficiently  far 
apart  to  be  considered  isolated.  First,  the  forces  acting  upon  several  structural  members  will 
be  determined.  The  total  force.  Fj(0,  S),will  be  a summation  of  the  drag  force,  Fp{0,  S), 
and  inertia  force,  F|(0,  S),  components  at  any  particular  phase  angle.  Each  component 
will  be  presented  graphically;  the  components  will  then  be  added  to  establish  the  total  force, 
and  the  maximum  force  acting  upon  each  member  will  be  obtained  from  the  graph. 

Forces  on  Member  "a  ” 

In  the  ca.se  of  the  outrigger.  Member  a,  the  drag  force  is  given  by: 


Cj^pD 

2 


S 

a 

u|u|dS  ' 

0 


where  IJ  is  tlie  piling  diameter,  Sg(=  20.5')  is  the  height  of  the  outrigger  above  bottom, 
and  p = mass  density  of  sea  water,  1.99  sluggs/foot’ . To  determine  Fp(0,  S^),  select  the 
tabulated  dimensionless  drag  value  for  the  force,  F^(0,S^),  at  depth  S|,/h  = 0.5  from 
Table  V and  multiply  the  dimensionless  force  by : 


CjjpD(H/T) 

2 


CjjPD(H/T) 

2 


1.05(1.99)  (6)  (31.78/20)  Ml 
2 


[648.9  lbs  ] 
[0.6489  kips] 
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Tin*  inertia  force  on  Vlemla'r  a is  {jiven  by: 


Fi(e,s^)  = 


CjjPttD' 


Du 

Dt 


dS'. 


To  determine  Fy(0,  S^),  stdeet  the  tabulated  value  of  the  dimensionless  inertia  force, 
Fj(0,  S),  for  a relative  depth  S^/h  = 0.5  from  Table  VI  and  multiply  the  dimensionless 
force  by: 


274.9  lbs 
0.2749  kips 

The  total  force  will  be  determined  by  summation  of  F|(0,  and  fpiO,  S^)  at  each  phase 
angle,  0.  The  force  calculations  are  summarized  in  Table  G and  the  forces  are  plotted  in 
Figure  31. 


TABLE  G 

Horizontal  Wave  Forces  on  Member  "a" 


ec) 

0 

10 

20 

30 

50 

75 

100 

130 

180 

36.31 

29.00 

14.60 

4.30 

- 0.04 

-1.14 

-1.54 

-1.62 

-1.60 

Fjjdcips) 

23.56 

18.81 

9.47 

2.79 

- 0.03 

-0.74 

-1.00 

-1.05 

-1.04 

'■l' 

0.0 

22.59 

36.36 

36.63 

17.25 

3.76 

0.67 

0.12 

0.0 

Fj  Ocips) 

0.0 

6.21 

10.00 

10.07 

4.74 

1.03 

0.18 

0.03 

0.0 

F^Uips) 

23.56 

25.02 

19.47 

12.86 

4.71 

0.29 

-0.82 

-1.02 

-1.04 

Forces  on  Member  “6” 

Next,  consider  the  horizontal  forces  acting  on  the  main  support  piling.  In  this  case,  the 
forces  are  integrated  from  0 to  h + t)(0).  To  determine  Fp(0),  multiply  the  tabulated 
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■e  31.  Horizontal  wave  forces  on  member 


value  for  the  dimensionless  total  drag  force,  F/>  (0)  (indicated  “Surface”  in  Table  V)  by 
the  same  constant  as  for  Member  a,  i.e., 


C pD(H/T)  "h 

2 “ 0.6849  kips 


Similarly,  Fj(0)  is  found  by  multiplying  the  tabulated  value,  F/(0),  indicated  “Surface”  in 
Table  VI  by; 

C^pttD^  (H/T^lh 

— 2 = 0.2749  kips 

The  calculated  forces  are  summarized  in  Table  H and  are  plotted  in  Figure  32. 


TABLE  H 

Horizontal  Wave  Forces  on  Member  "b" 


ec) 

0 

10 

20 

30 

50 

75 

100 

130 

180 

242.39 

119.80 

37.00 

7.72 

- 0.25 

-2.19 

-2.84 

-2.95 

-2.92 

F^(lcips) 

157.3 

77.7 

24.0 

5.0 

- 0.2 

-1.4 

-1.8 

-1.9 

-1.9 

0.0 

112.13 

113.47 

84.55 

30.12 

6.08 

1.03 

0.27 

Pj(kips) 

H 

30.8 

31.2 

23.2 

8.3 

■ 

0.3 

0.1 

mm 

P,j.(kip8) 

157.3 

108.  5 

55.2 

28.2 

8.1 

1 

B 

B 

Forces  on  Member  “c" 

Finally,  consider  structural  Member  c,  the  fender.  The  computation  for  this  member  is  a 
combination  of  the  two  previous  methods  since  it  is  sometimes  over-topfied  by  the  wave. 
The  forces  are  integrated  from  Sp^  = 32.8  ft  to  Sp^  = 45.1  ft;  therefore,  the  force  acting  on 


Figure  32.  Horizontal  wave  forces  on  menriber 


an  iniapnary  piliiif;  up  to  the  bottom  of  the  fender  is  subtracted  from  a similar  term  for  the 
top  of  the  fender.  The  dimen.sionless  forces  are  olitained  by  .subtracting  the  dimen.sionless 
force  components  |H'rtaining  to  the  bottom  of  the  member  from  those  pertaining  to  the  top. 
If  the  top  of  the  member  is  submerged,  the  value  at  = 1 .1  .should  be  used;  for  times  that 
the  top  is  not  .subm<*rged,  the  value  indicated  “.Surface”  should  be  employi^d  for  Note 
that  the  selection  of  the  proper  value  for  the  member  upper  elevation  follows  readily  from 
the  tables:  the  values  at  S|,^  = 1.1  are  used  at  phase  angles  where  they  are  tabulated 
(0<0<20°)  and  the  values  labeled  “Surface”  arc  u.s<*d  for  the  remaining  pha.se  angles 
<0  < 1R0°). 

Summari/.ing,  for  each  phase  angle,  the  net  dimensionless  force  components  on  Member  c 
are  obtained  by: 


where  the  subscripts,  N.  U and  L indicate  net,  upper  and  lower.  The  dimensionalizing 
constant  for  drag  force  for  the  member  is  calculated  (recalling  that  D = 3') 

C-pD(H/T) 2 h 

-2 — 2 = 0.3  245  kips 

and  for  the  inertia  force  component 

C„PttD^  (H/T2)h 

^ = 0.0687  kips 

The  required  calculations  arc  summarized  in  Table  I and  the  results  are  shown  in  Figure  33. 

The  maximum  horizontal  wave -induced  forces  are  now  available  for  the  design  wave,  and 
may  be  used  in  further  design  analysis.  They  are  summarized  in  Table  J. 

Moments  on  Member  “a” 

The  moments  due  to  the  wave  forces  acting  on  the  structure  are  also  essential  in  design. 
For  any  member,  the  moment  about  the  mudline  is  defined  as; 


fS2 

■ Sj 

fS2 

S dF  (0,s)  = 

S dF  (0,S)  + 

S dFj(0,S) 

Si 

Si 

^Si 

= Mp(0)  + Mj(0) 


65 


TABLE  I 

Horizontal  Wave  Forces  on  Member  "c” 


0 

10 

20 

30 

50 

75 

100 

130 

150 

99.73 

75.87 

33.17 

7.72 

- 0.25 

-2.19 

-2.84 

-2.95 

-2.92 

63.34 

49.68 

23.72 

6.40 

- 0.14 

-1.87 

-2.48 

-2.59 

-2.56 

36.39 

26.19 

9.45 

1.32 

- 0.11 

-0.32 

-0.36 

-0.36 

-0.36 

11.81 

8.50 

3.07 

0.43 

- 0.04 

-0.10 

-0.12 

-0.12 

-0.12 

0.0  65.78  96.88  84.55  30.12  6.08  1.03  0.27 


0.0  40.55  62.97  60.49  26.23  5.53  0.96  0.22 


PjOcips) 


25.23  33.91  24.06  3.89  0.55  0.07  0.05 


1.73  2.33  1.65 


F^(Itips)  11.81  10.23  5.40  2.08 


0.27  0.04 


0.23  -0.06  -0.11  -0.11  -0.12 


where 


and 


Md(0) 


PCdD 

~5“ 


>82 

Si 


s 


u I u I dS 


Mj(0,S)  = I ^ S ^ dS 


TABLE  J 

Summary  of  Maximum  Wave  Forces  on 
Several  Platform  Components 


Phase  Angle  of 
Maximum  Force , 

F (kips) 

Member 

0m 

n\ 

max 

a 

70 

25.1 

b 

1“ 

160 

c 

1° 



12.3 

Note:  Phase  angles  and  maximum  forces  obtained  by  interpolation 
from  Figures  31 , 32  and  33. 


Consider  the  total  moment  about  the  mudline  on  the  outrigger  (Member  a).  In  this  case 
S,  = 0,  and  Sj  = = 0.5  h.  To  determine  the  drag  moment,  Mp(0),  multiply  the 

dimensionless  tabulated  value  for  the  drag  moment,  M^(0),  listed  at  depth  S^/h  = 0.5  in 
Table  VH.by: 


CjjPD(H/T) 

2 


26.606  for  Mq  in  ft-lbs 

26.606  for  in  ft-kips 


Similarly,  multiply  M|(0)  listed  at  depth  S^/h  = 0.5  in  Table  VIII  by: 


C^ptrD*  (H/T^)h* 

3 


11.272  for  Mj  in  ft-lbs 

11.272  for  Mj  in  ft-kips 
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to  obtain  M|(0).  Thew  moments  arc  added  to  obtain  My(0),  as  siiown  in  Table  K and 
Figure  34. 


TABLE  K 

Wave  Moments  (About  Mudline)  on  Member  "a" 


ec) 

0 

10 

20 

30 

50 

75 

100  130 

180 

“d 

9.31 

7.40 

3.67 

1.05 

- 0.01 

- 0.29 

- 0.39  - 0.40 

- 0.40 

Motft- 

kips) 

247.7 

196.9 

97.6 

27.9 

- 0.3 

- 7.7 

-10.4  -10.6 

-10.6 

0.0 

5.85 

9.32 

9.26 

4.25 

0.92 

0.16  0.03 

0.0 

kips) 

0.0 

65.9 

105.1 

104.4 

47.9 

10.4 

1.8  0.3 

0.0 

kips) 

247.7 

262.8 

202.7 

132.3 

47.6 

+ 2.7 

- 8.6  -10.3 

-10.6 

Moments  on  Member  “6” 

Next  consider  the  moment  on  tbe  main  structural  piling  (Member  b).  The  limits  of 
integration  are  from  0 to  h + r)(d).  Therefore,  take  the  tabulated  values  labeled  “Surface” 
from  Table  VII,  [M;,(0)],  and  Table  VIII,  [Mj(0)],  and  multiply  by: 

Cj-pD(H/T) 

2 = 26.606  for  Mp  in  ft-kips 

and 

Cj-PTTD^  (H/T*)h2 

^ = 11.272  for  Mj  in  ft-kips 

in  order  to  obtain  and  My(0).  The  two  moments  are  added  to  obtain  Mj(0)  as 

indicated  in  Table  L and  plotted  in  Figure  35. 

Moments  on  Member  "c" 

The  fender  has  the  same  limits  of  integration  for  moment  calculation  as  for  the  force 
calculation  and  is  determined  in  a similar  manner.  However,  the  tablulatcd  moments, 
M^(0,  .S),  and  My(0,S),  are  taken  from  Tables  VII  and  VIII.  The  total  moment  acting  on 
the  fender  is  found  by:  Mj'(0)  = Mp(0)  + M^(0).  The  calculations  are  summarized  in 
Table  M and  are  plotted  in  Figure  36. 
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TABLE  L 


Wave  Moments  (About  Mudline)  on  Member  "b" 


ec) 

0 

10 

20 

30 

50 

75 

100 

130 

180 

268.1 

102.6 

23.0 

3.6 

- 0.2 

- 1.0 

- 1.3 

- 1.3 

- 1.3 

M^(ft- 

klps) 

7133 

2730 

612 

96 

- 5 

-27 

-35 

-35 

-35 

0.0 

101.7 

78.5 

47.5 

13.5 

2.5 

0.4 

0.1 

0.0 

kips) 

0.0 

1146 

885 

535 

152 

28 

5 

1 

0.0 

kips) 

7133 

3876 

1497 

631 

147 

1 

-30 

-34 

-35 

TABLE  M 

Wave  Moments  (About  Mudline)  on  Member  "c" 


ec) 

0 

10 

20 

30 

50 

75 

100 

130 

180 

0 

61.94 

46.01 

18.59 

3.63 

- 0.18 

-1.04 

-1.31 

-1.35 

-1.33 

n 

27.04 

20.94 

9.61 

2.40 

- 0.08 

-0.77 

-1.00 

-1.04 

-1.02 

131 

34.90 

25.07 

8.98 

1.23 

- 0.10 

-0.27 

-0.31 

-0.31 

-0.31 

(ft-kipa) 

464 

334 

119 

16 

- 1 

-4 

-4 

-4 

-4 

0.0 

41.87 

59.20 

47.47 

13.45 

2.52 

0.40 

0.14 

0.0 

H 

0.0 

17.66 

26.76 

24.82 

10.04 

2.05 

0.34 

0.10 

0.0 

n 

in 

24.21 

32.44 

22.65 

3.41 

0.47 

0.06 

0.04 

(ft-klps) 

68 

91 

64 

10 

1 

0 

0 

0.0 

(ft-kipa) 

464 

402 

210 

80 

9 

-3 

-4 

-4 

-4 

The  maximum  calculated  forces  and  moments  on  the  three  platform  members  due  to  the 
design  wave  are  summarized  in  Table  N. 


TABLE  N 

Summary  of  Maximum  Wave  Forces  and  Moments 


Member 

0 

e 

Mt(0,S) 

(ft-kips) 

a 

70 

25.1 

B 

267 

b 

1“ 

160 

714  0 

c 

1“ 

12.3 

B 

475 

Example  2— Wave  Characteristics,  Kinematics  and  Pressure  Fields 

TTiis  example  describes  the  use  of  the  tables  for  calculating  various  parameters  associated 
with  a periodic  wave.  These  parameters  include  the  wavelength  and  the  kinematic  and 
pressure  fields. 

Wavelength 

The  wavelength  is  presented  in  dimensionless  form  in  Table  XI  of  the  sample  output  and 
is  determined  as  follows: 


L 


L ' 


For  example,  for  the  same  wave  considered  in  Example  1,  L'*=  0.422  and  T - 20  seconds. 
TTie  wavelength  is  therefore: 

L » 5.12  (20)^(0.422)  = 864.3  ft 


Wave  Profile 

The  dimensionless  wave  profile,  ri\d),  is  tabulated  in  each  of  Tables  I through  IX  and  is 
defined  as: 


n ' (0) 


n(9) 

-n — 
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therefore 


n (0)  - n ' (0)  • H 

Hw  wave  profUe  calculatioii  for  Caae  4-D  ia  niminarized  in  TaMe  O and  is  plotted  in 
FIfure  37.  Note  that  q is  an  even  function  of  0. 

WaierHrtid*  Kmemmtie$ 

Tlie  water  particle  kinematics  will  be  calculated  for  Case  4-D  as  presented  in  the  sample 
output.  These  kineinatics  will  be  calculated  for  mid-depth  (ix.,  20.5  feet  above  the  bottom). 
Hie  dimensionless  forms  of  these  variables  are  presented  in  Tables  I through  IV  of  the 
ma^de  output.  The  dimensionlem  water  particle  velodties  are  defined  as: 

u'(0,S)=H<0^ 


and 


...  » /a  e»  = 'l^(0  »S) 


and  the  dimenaionlem  water  particle  total  accelerations  are  defined  as: 


and 


Du 

- 

■ 


Dw 
- Dt 
“ BTTT 


Note  that  these  are  functions  of  $ and  S,  however,  for  convenksnce,  the  dependence  has 
not  been  indicated  in  the  above  expressions.  The  calculations  of  the  water  particle  velocities 
and  accelerations  over  the  ranfe  0*  < 9 < 180*,  are  also  summarixed  in  Table  O and  plotted 
in  FIfure  37. 

It  will  be  noted  that  in  the  tables  of  wave  functioiu,  the  variables  are  only  presented  for 
phase  angles  ranging  between  0*  and  180*.  All  of  the  variables  sre  either  symmetrical  or 
antiaymmetrical  about  a phase  angle  of  0*.  The  variables  that  are  symmetrical  include:  the 
water  surface  profile,  the  horixontal  component  of  water  particle  velocity  and  the  vertical 
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TABLE  O 


* 

\ 


component  of  water  particle  acceleration.  The  antisymmetrical  variables  include  the  vertical 
component  of  velocity  and  the  horizontal  component  of  velocity  and  the  horizontal 
component  of  water  particle  acceleration. 

Dynamic  Pressure 

The  dynamic  pressure  also  was  calculated  at  a distance  of  20.5  feet  above  the  bottom. 
The  dimensionless  form  of  this  variable  is: 


IH 

2 


and  is  presented  in  Table  IX  of  the  sample  output.  The  calculations  are  summarized  in 
Table  0 of  this  report  and  presented  in  graphical  form  in  Figure  37.  Note  that  p^  is  an  even 
function  of  0. 

Example  3— Free  Surface  Boundary  Condition  Errors 

The  free  surface  boundary  condition  errors  and  the  reason  for  examining  and  tabulating 
these  errors  have  been  described  in  Section  II.  To  illustrate  the  use  of  tables  to  calculate  the 
free  surface  boundary  condition  errors,  both  the  distributed  errors  on  the  free  surface  and 
the  root-mean-square  and  maximum  errors  as  gross  measures  of  these  errors  will  be 
presented.  The  distributed  kinematic  and  dynamic  free  surface  boundary  condition  errors 
are  presented  in  Table  X,  Items  1 through  4 of  the  sample  output  and  the  root-mean-square 
errors  and  maximum  errors  are  presented  in  Table  XI,  Items  10  through  13. 

Distributed  Boundary  Condition  Errors 

The  calculations  of  the  distributed  boundary  condition  errors  are  presented  in  Table  P 
and  Figure  38.  The  kinematic  free  surface  boundary  condition  errors  as  defined  and 
presented  in  the  wave  tables  (Table  X,  Items  I and  2)  are  in  dimensionless  form.  However, 
the  dynamic  free  surface  boundary  condition  errors  (Table  X,  Items  3 and  4 of  wave  tables) 
are  dimensional  as  illustrated  in  the  sample  calculations  accompanying  Table  P.  The 
calculations  of  the  root-mean-square  (RMS)  and  maximum  kinematic  dynamic  free  surface 
boundary  condition  errors  are  presented  below. 

Overall  Kinematic  Free  Surface  Boundary  Condition  Errors 

The  RMS  kinematic  free  surface  boundary  condition  errors  are  presented  as  Item  10  in 
Table  XI,  i.e., 

0.0475  (Linear  Wave  Theory) 

0 . 0 (Stream-function  Wave  Theory) 
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Free  Surface  Boundary  Condition  Errors 


e2(ft)  0.918  -0.356  -0.343  -0.124  0.022  0.006  0.064  0.041  0.010 


The  maximum  KFSBC.  error  is  obtained  from  Item  12  of  Table  XI, 

Iciljng^x  “ 0.0856  (Linear  Wave  Theory) 

1 e 1 1 =0.0  (Stream-function  Wave  Theory) 

max  ^ 

Owrall  Dynamic  Free  Surface  Boundary  Condition  Errors 

The  RMS  DFSIU;  errors  are  presented  in  dimensionless  form  as  Item  11  in  Table  Xl,i.e., 


/ 


0.0241 


= 0.765  ft  J 


(Linear  Wave  Theory) 


0.0048 


(Stream-function  Wave  Theory) 


= 0.153  ft 

The  maximum  DFSBt;  errors,  obtained  from  Table  XI,  Item  13  are: 


£2 


-niax_  = 0.0385 


lez I = 1.224  ft  J 

max 

= 0.0289 

H 

|e2|  = 0.918  ft 

max 


(Linear  Wave  Theory) 


(Stream -function  Wave  Theory) 


In  the  interpretation  of  the  boundary  condition  errors  in  accordance  with  the  discussion  in 
Section  II,  if  the  Iwundary  condition  errors  for  any  given  tlieory  were  found  to  be  generally 
better  than  for  the  Stream-function  theory,  then  it  could  be  concluded  that  at  least  the 
analytical  validity  of  that  wave  theory  would  be  better,  and  (as  discussed  earher)  there  is 
evidence  that  the  analytical  wave  theory  is  a good  indicator  of  the  experimental  validity  (or 
of  the  wave  phenomenon  in  nature). 
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F.^ampl#  4— Cakulation  of  Energy,  Momentum,  and  Energy  and  Momentum  Flnxea 

TTie  tabulations  of  average  potential,  kinetic,  and  total  energy  and  energy  fluxes  and 
average  momentum  and  momentum  fluxes  are  presented  in  Table  XI,  TTie  calculation  of 
these  quantities  in  dimensional  form  is  straightforward  and  wiD  simply  be  presented  without 
discussion. 

Average  Potential  Energy  (Table  XI,  Item  2) 


PE  = 0.213(8080)  = 1721  ft-lb/ft^ 


Average  Kinetic  Energy  (Table  XI,  Item  3) 


TyHVST 


= 0.254 


KE  = 0.254(8080)  = 2052  £t-lb/ft^ 


Total  Energy  (Toble  XI,  Item  4> 


" (yttVS)  “ 

TE  = 0.467(8080)  * 3773  ft-lb/ft* 


Energy  Flux  (Table  XI,  Item  5) 

F'  = ■ . = 0.447 

TE  yH*  LI 

I 

F-,„  = 0.447(349166)  - 156077  £t-lbs/ (ft-sec) 

TE 


Group  Velocity  (Table  XI,  Item  6} 


“ (L^T)  “ 


0.957(43.21)  « 41.36  £t/8ec 
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Average  Momentum  (Table  XI,  Item  7) 


M' 


0.505 


M = 

Average  Momentum  Flux 


F'  = 
m 


0.505(187)  * 94.42  Ib-sec/ft* 

in  IFave  Direction  (Table  XI,  Item  8) 


0.603 

IVJ 


0.603(8080)  = 4872  Ib/ft 


The  average  momentum  flux  has  been  recognized  in  recent  years  as  an  important  dynamic 
quantity  and  is  related  to  wave  setup  within  the  surf  zone  and  also  is  an  important  factor  in 
the  longshore  transport  of  littoral  material. 

Average  Momentum  Flux  Transverse  to  IFave  Direction  (Table  XI,  Item  9) 


0.156 


= 0.156(8080)  = 1260 


From  the  momentum  flux  components  presented,  it  is  possible  to  obtain  any  component 
of  the  radiation  stress  tensor  (Bowen,  1969). 

Example  5— Free  Surface  Breaking  Parameters 

The  free  surface  breaking  parameters  as  defined  by  Equations  (48)  and  (49)  are  based  on 
two  stability  considerations.  The  kinematic  free  surface  breaking  parameter  is  defined  in 
terms  of  the  speed  of  a water  particle  on  the  surface  at  the  crest  relative  to  the  wave  form 
speed.  If  this  parameter  should  equal  unity,  then  the  wave  is  regarded  as  unstable  due  to 
kinematic  considerations.  The  dynamic  free  surface  breaking  parameter  is  defined  as  the 
ratio  of  the  vertical  acceleration  of  a water  particle  on  the  surface  at  the  wave  crest  relative 
to  the  acceleration  of  gravity.  The  interpretation  is  that  if  this  parameter  should  equal  unity , 
then  the  pressure  immediately  under  the  crest  would  be  zero  and  if  the  parameter  should 
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exceed  unity,  tlicn  according  to  the  equations  of  motion,  the  pressure  beneath  the  wave 
crest  would  be  negative  which  is  unrealistic  and  would  indicate  an  unstable  water  surface. 

It  should  be  noted  that  the  theory  employed  in  the  study  is  composed  of  a finite  series  of 


terms.  To  adequately  define  an  instability  formally,  it  may  be  necessary  to  extend  the 
representation  to  include  an  infinite  number  of  terms.  The  results  presented  here  for  the 
free  surface  breaking  parameters  should  be  interpreted  accordin^y.  For  the  sample  output 
(Case  4-D,  Table  XI,  Item  14)  shows  that  the  kinematic  free  surface  breaking  parameters  for 
the  linear  and  Stream -function  representations  are  0.429  and  0.733,  respectively.  The 
corresponding  values  (Table  XI,  Item  15)  for  the  dynamic  free  surface  breaking  parameter 
are  0.0409  and  0.286,  respectively.  The  wave  height  associated  with  this  case  is 
approximately  0.78  of  the  depth  and  according  to  the  McGowan  criterion,  the  wave  would 
be  breaking. 

Example  6— Combined  Shoaling-Refraction 

The  shoabng-refraction  results  were  not  tabulated,  but  are  presented  for  various 
deepwater  directions  in  ^aphical  form  as  Figures  25  through  29  of  this  report. 

Example  6-a 

Consider  a deepwater  wave  propagating  over  bathymetry  characterized  by  straight  and 
parallel  contours;  the  deepwater  wave  conditions  considered  are: 

Ho  = 11.52  ft 
T =15  sec 
oo  = 40® 

Suppose  that  we  wisli  to  find  the  wave  height  and  direction  in  a water  depth  of  30  feet  and 
also  the  wave  height,  water  depth  and  wave  direction  at  breaking.  Figure  28  is  applicable  for 
a deepwater  wave  direction  of  40**.  The  deepwater  wavelength  is  calculated  as: 

therefore 

Hs 


TTie  line  for  = 0.01  is  simply  followed  to  the  left  to  the  intersection  with  h/L^,  = 

0.0260.  At  this  intersection, 

= 0.0119 

liO 

H - (0.0119) (1152)  = 13.71  ft 
a » 17" 

The  second  part  of  the  example  requires  the  breaking  depth,  height  and  angle.  For  this,  the 
H^/L^j  = 0.01  curve  intersects  the  breaking  curve  at: 

= 0.0190 
Lo 

therefore 

“b 

~ » 0.0147 
Lo 

“b  “ 

therefore 

Hg  = 0.0147(1152)  - 16.9  ft 
hg  = 0.0190(1152)  » 21.9  ft 

Example  6-b 

Suppose  that  a wave  is  observed  in  transitional  depths  and  it  is  desired  to  determine  the 
height  at  deep  water,  breaking,  or  any  depth  of  interest.  For  this  example,  the  values 
of  H/L^  and  h/L^  are  calculated  from  the  observed  wave  height  and  period  and  water 
depth.  If  the  observed  direction  corresponds  to  one  of  the  graphs  available,  then  one 
proceeds  as  before  in  Example  6-a.  If  the  observed  point  is  not  in  accordance  with  any  of 
the  graphs  available,  then  an  interpolative  procedure  is  required.  As  an  example,  consider 
the  following  observed  wave  characteristics 

H « 20  ft 
h " 60  ft 
T » 12  sec 
o « 11" 
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and  it  is  desired  to  ealculate  the  wave  height  and  direction  in  a water  depth  of  40  feet.  From 
the  observed  information 


Lo  = 737.3 

H/Lo  = 0.0271 

h/Lo  = 0.0814  (h  = 60  ft) 

h/Lo  = 0.0542  (h  = 40  ft) 


Examining  the  available  figures,  it  is  seen  that  the  deepwater  wave  direction  is  between 
10°  and  20°.  As  a close  approximation,  the  problem  is  solved  for  = 10°  and  = 20°, 
and  the  desired  results  obtained  by  interpolation.  For  = 10°,  from  Figure  26,  a line 
passing  through  H/L^  =0.0271,  h/L^=  0.0814  is  sketched  with  the  same  approximate 
shape  as  those  for  = 0.02  and  0.04  to  determine  H/^  = 0.033  and  a = 6.2°  for 

h/Lj,  = 0.0542.  The  corresponding  values  for  = 20°  are  H/L^  = 0.031  and  a = 12°.  The 
procedure  is  shown  graphically  in  Figure  39  for  = 10°.  Because  ft  = 10°  and  20°,  the 
O values  corresponding  to  h/L^  = 0.0814  and  H/L^^  = 0.0271  are  6.8°  and  13°  respectively, 
and  the  desired  a for  these  conditions  is  11°,  the  values  of  H/L^  and  a for  h = 40  feet  may 
be  determined  by  linear  interpolation  as: 

= 0.033  + i'lsi  I 5:5»y  = 0*032 

or 

H = (737.3) (0.032)  = 23.6  ft 
and 

a = 6.2»  + I l\ll\  (11°  - 6.8“)  = 10.1» 

Dissipative  mechanisms  such  as  percolation  and  bottom  friction  are  not  included  in  these 
results,  and  in  many  cases  these  mechanisms  will  be  of  greater  significance  than  the 
nonlinear  effects  on  the  celerity  and  group  velocity  which  represent  the  difference  between 
the  results  presented  here  and  the  linear  wave  theory. 

Example  7— Use  of  Tables  for  Nontabulated  Wave  Conditions 

Most  of  the  previous  examples  have  been  presented  for  wave  conditions  which  were 
available  as  one  of  the  40  tabulated  cases,  i.e.,  Case  4-D.  It  is  anticipated  that  the  tabulations 
will  be  used  primarily  for  preUminary  design,  and  therefore  that  the  40  cases  may  provide 
adequate  information  for  this  purpose  without  interpolation.  Final  design  of,  for  example,  a 
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H/L„  = 0.0271  to  h/L„  = 0.0542 


platform  supported  by  battered  piling  would  probably  be  carried  out  by  establishing  a 
Stream-function  (»r  other  wave  theory  representation  for  the  particular  wave  conditions 
selected  for  design. 

On  occasion,  it  may  be  desired  to  interpolate  between  the  ca8f;8  presented  in  the  tables 
for  wave  conditions  that  are  substantially  different  from  any  of  the  40  cases.  Several 
numerical  and  graphical  interpolation  methods  were  explored  with  a goal  of  obtaining  a 
simple  method  which  yielded  reasonably  accurate  results.  Because  most  wave  variables  of 
interest  are  nonlinear,  numerical  schemes  which  used  Unear  interpolation  proved  to  be 
inaccurate.  The  best  procedure  was  found  to  be  a simple  graphical  procedure  which 
generally  yields  results  within  5 percent. 

Method 

The  method  uses  the  tabulated  parameters  of  interest  for  the  H/Hg  values  above  and 
below  the  value  of  interest  at  the  two  lower  and  two  higher  h/L^  tabulated  values;  in  all 
for  each  parameter  desired,  the  interpolated  value  is  based  on  values  of  that  parameter  for 
eight  tabulated  wave  conditions.  The  method  is  outUned  in  the  following  paragraphs  and 
illustrated  by  two  examples. 

Suppose  that  the  wave  height,  period,  and  water  depth  selected  for  design  are  H^,  Td- 
and  h^.  The  design  wave  steepness  and  relative  depth  are  calculated  as: 


Wave  Steepness; 

Lod 


Relative  Depth:  = 

Ijod 


where 


Loi 


= -2- 
2it 


The  relative  depth  and  wave  steepness  are  plotted  on  Figure  40  to  establish  which  wave 
cases  should  be  used  for  design.  For  the  example  shown,  = 0.086  and 

h/Lg^  = 0.313.  This  point  falls  between  H/Hg  values  denoted  as  B and  C (i.e.,  50  and  75 
percent  of  breaking  heights,  respectively)  and  between  tabulated  h/L^  values  denoted  as 
Cases  7 and  8.  The  interpolation  would  therefore  be  based  on  the  tabulated  parameter  of 
interest  for  Cases  6-B,  6-C,  7-B,  7-C,  8-B,  8-C,  9-B,  and  9-C. 

The  interpolation  proceeds  as  follows.  An  auxiliary  plot  is  made  of  the  variable  of 
interest,  e.g.,  the  total  dimensionless  drag  force  at  0 = 0°  [denoted  Fp(0‘*,  Surf.)).  This  plot 
provides  a continuous  distribution  of  Fj)(0°,  Surf.)  versus  h/L^  for  relative  breaking 
heights  B and  C.  Interpolated  F^  values  are  then  obtained  from  the  auxiliary  plot  for 
the  h/Lg  design  value  (0.313).  The  interpolation  for  the  design  wave  steepness  requires 
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Figure  40.  Interpolation  aid 
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measuring  (Figure  40)  the  vertical  linear  distance  from  the  B and  C lines  to  the 
design  H/L^  of  interest;  denote  these  values,  A,  and  respectively.  Weighting 
factors,  W,  are  then  established  as: 


W, 


Aa 


Ai  + A2 


W. 


O. 


'U  Ai  + Aj 
The  interpolated  value  is  finally  determined  as: 


(50) 


(Fd)d  = Wl  (FiJ)j^  + Wu  (Fi;)„ 


where  the  subscripts,  D,  L and  U outside  the  parentheses  denote:  “Design,”  “Lower” 
(Case  B),  and  “Upper”  (Case  C). 

Example  7-a— Numerical  Illustration  of  Interpolation  Procedure 
Consider  the  following  wave  conditions  selected  for  design 


which  yield 


Hp  * 44  ft 
Tjj  » 10  sec 

hjj  » 160  ft 


^ T*  » 512  ft 


0.313 

0.0859 


and  suppose  that  we  require  the  maximum  dimensionless  drag  force  on  a piling  that  extends 
from  the  bottom  to  above  the  crest  level.  This  maximum  value  would  occur  at  0 = 0°  and  is 
the  value  labeled  “SURFACE”  in  the  tabulations.  Plotting  of  the  wave  steepness  and  relative 
depth  on  Figure  40  indicates  that  the  design  values  are  spanned  by  Cases  7-B,  7-C,8-B  and 
8-C.  In  accordance  with  the  preceding  section  the  values  of  Fp(0“,  Surf.)  for  Cases  6-B,  6-C, 
7-B,  7-C,  8-B,  8-C,  9-B,  and  9-C  are  required  for  interpolation  and  are  summarized  in 
Table  Q. 
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TABLE  Q 


Summary  of  Fj^(0®,  Surf. ) Required  for  Exeunple  7-a 


Case 

FpIO®,  Surf.) 

6-B 

22.37 

6-C 

28.79 

7-B 

8.60 

7-C 

11.31 

8-B 

2.71 

8-C 

3.53 

9-B 

1.33 

9-C 

1.72 

The  values  in  Table  Q are  presented  as  an  auxiliary  plot  in  Figure  41.  Interpolation  at  the 
design  h/L^  of  0.313  yields  the  following  values  of  F^  for  relative  breaking  of  50  and  75 
percent  respectively. 


Relative  Breaking  of  50  percent  (Line  B):  (F^))^  = 4.90 
Relative  Breaking  of  75  percent  (Line  C):  (P'D)fj  - 6.10 

To  interpolate  to  the  design  H/L^,  the  distances  A,  and  A,  are  measured  from 
Figure  40.  For  this  example,  these  are  found  to  be: 


Li  • 0.11  in 
62  » 0.32  in 

Tlie  weighting  values  are  then  (Equation  50) 

w _ 62 

L SI  + A2 

u K A 1 

U Sx  + A2 


0.744 


« 0.256 


and  the  inti'rpolatrd  value  of  F},  is: 


= Wl(Fi^)j^  + Wu(Fi5)y 

= (0.744)  (4.90)  + (0.256)  (6.10) 

= 5.21 

To  evaluate  this  interpolated  value,  a Stream-function  solution  was  developed  for  the 
conditions  of  interest,  and  from  the  actual  solution  was  found  to  be  5.04  or  a difference 
of  about  ,3.4  jiereent. 

More  comprehensive  evaluations  of  the  accuracy  of  the  interpolation  method  are 
presented  in  the  next  example. 

Example  7-b— Assessment  of  the  Interpolation  Method 

To  present  a more  extensive  evaluation  of  the  accuracy  of  the  interpolation  method,  two 
sf)ecial  eases  (one  shallow-water  and  one  deepwater)  were  .selected  for  evaluation.  The  wave 
characteristics  for  these  two  cases  are  presented  in  Table  R. 


TABLE  R 

Wave  Characteristics  Selected  for  Accuracy 
Evaluation  of  Interpolation  Method 


Case 

Wave  Height, 
H(ft) 

Wave  Period, 

T (sec) 

Water  Depth, 
h(ft) 

S-1 

(Shallow  Water) 

19 

20 

30 

S-2 

(Deep  Water) 

44 

10 

160 

Using  the  procedure  de8<Tibed,  interpolated  values  of  a number  of  quantities  of  interest 
were  develo|)cd  and  compared  with  values  obtained  by  Stream-function  solutions  at  the 
wave  conditions  of  interest.  Table  S presents  a summary  of  the  percentage  differences 
between  the  solution  and  interpolated  values. 

As  an  overall  statement  regarding  the  interpolation,  it  is  noted  that  Table  S indicates  that 
the  procedure  presented  generally  provides  results  which  are  within  5 percent;  however, 
differences  up  to  10  fiercent  could  occur.  One  final  comment  concerning  the  consistency  of 
the  tabulated  values  is  in  order.  In  preparing  the  auxiliary  plots,  it  was  usually  found  that  a 
bne  could  be  dr.'> wn  through  the  four  points  within  2 to  3 percent,  except  for  the  breaking 
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TABLE  S 


Svunmary  of  Percentage  Differences  Between  Values 
Determined  by  Streeun  Function  Solutions 
and  by  Interpolation 


Percentage  £ 

)if  ference^ 

Dimensionless  Varledjle^ 

Case  S-1 

Case  S-2 

u'(0*,  0.5);  Horiz.  Vel. 
Comp.,  Zero  Phase  Angle, 
Mid-Depth 

+3.9 

<1 

Fq(0°,  Surf),  Max.  Drag 
Force  Component,  Acting 
Over  Entire  Depth 

+6.7 

+3.4 

Fj(10®,  Surf),  Inertia 

Force  Component 

+1.3 

Not  Evaluated 

FjCTS",  Surf),  Inertia 
Force  Component 

Not  Evaluated 

-3.9 

Mi5(0**,  Surf),  Max.  Drag 
Moment  Component  About 
Mudline 

+4.5 

+3.6 

Mx(10^,  Surf),  Inertia 

Moment  Component 

+2.2 

Not  Evaluated 

Mi  (75®,  Surf),  Inertia 

Moment  Component 

Not  Evaluated 

-3.7 

Pd(0®,  0.5),  Dynamic 
Pressure  Component, 

Zero  Phase  Angle, 

Mid-Depth 

<1 

-2.4 

P]^(180®,  0.5),  Dynamic 
Pressure  Component, 

Trough  Phase  Position, 
Mid-Depth 

<1 

-2.8 

L',  Have  Length 

1.1 

<1 

TE',  Total  Energy 

-4.6 

-3.7 

TABLE  S — Continued 


Percentage  Difference!} 

Dimensionless  Variedale'^ 

Case  S-1 

Case  S-2 

Frpg,  Total  Energy  Flux 

-4.2 

+ 3.5 

M ' , Momentum 

-4.1 

-2.2 

Momentum  Flux  in 

Wave  Direction 

-3.7 

-2.6 

Fm  , Momentum  Flux 
^ Transverse  to 

Wave  Direction 

-1.7 

<1 

KFSBP,  Kinematic  Free 

Surface  Breaking 

Pareuneter 

• 

00 

+4.4 

DFSBP,  Dynsunic  Free 

Surface  Breaking 

Pareuneter 

1.4 

<1 

“Refer  to  Tables  D,  E,  and  F for  a more  complete  description  of  the  dimensionless  variables. 


bn  . n vM  _ Interpolated  Value  — Stream-function  Solution 

Percentage  Difference  = ^ 

Stream-function  Solution 


wave  height,  H/Hg  = l.O  in  which  case  the  maximum  deviations  could  amount  to 
±5  percent.  Tlie  probable  explanation  f«>r  this  deviation  is  that:  (l)the  calculated  wave 
heiglits  for  the  tabulated  cases  were  allowed  to  deviate  from  the  desired  values  by  1 |)ercent, 
and  (2)  the  different  orders  lo  represent  different  cases  could  cause  a difference  in 
kinematics  of  I to  2 percent.  The  effects  noted  above  could  conceivably  amount  to 
deviatiuns  of  ±5  percent  for  those  variables  which  are  inherently  nonlinear,  e.g,,  drag  forces 
or  wave  breaking  parameters. 

This  completes  the  section  illustrating  the  use  of  the  wave  tables.  It  should  be  recognized, 
however,  that  only  the  more  simple  examples  have  been  presented  and  that  the  tables  can  be 
effectively  applied  to  the  solution  of  situations  which  are  considerably  broader  and  more 
complex  than  those  examined  in  this  section. 

VI.  SUMMARY 

This  report  presents  the  results  of  an  investigation  which  has  demonstrated  that  the 
Stream-function  wave  theory  provides  a generally  better  representation  of  periodic  wave 
phenomena  than  other  wave  theories  examined.  As  a result  of  this  indication,  tables  have 
been  prepared,  based  on  the  Stream-function  wave  theory,  that  include  parameters  which 
should  be  an  aid  in  preliminary  offshore  design.  The  tables  also  include  parameters  which 
are  presently  of  greatest  interest  to  researchers. 

Because  of  its  simplicity,  the  linear  wave  theory  is  widely  used  for  many  calculations  over 
all  ranges  of  relative  depth.  This  study  has  identified  that,  for  a number  of  variables,  there 
are  substantial  differences  between  the  linear  and  Stream -function  wave  theories.  Althougli 
this  point  has  not  been  amplified  in  this  report,  inspection  of  the  tables  will  substantiate  this 
conclusion.  The  identification  of  these  differences  should  be  of  assistance  in  planning 
experimental  programs  to  provide  definitive  research  results. 

If  the  set  of  tables  is  extensively  appUed,  as  is  hoped,  undoubtedly  the  users  will  note 
shortcomings,  omissions  or  develop  recommendations  directed  toward  the  improved 
usefulness,  applicability,  or  efficiency  of  the  tables.  The  author  would  welcome  information 
of  this  type  so  that  future  work  may  benefit  by  as  wide  a range  of  user’s  needs  as  possible. 
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APPENDIX  1 


NUMERICAL  SOLUTION  OF  STREAM  FUNCTION  PARAMETERS 
Introduction 

Thia  appendix  outlines  the  method  of  determining  numerical  values  for  the  parameters  in 
the  general  form  of  the  Stream-function  solution.  The  numerical  solution  requires  the  use  of 
a reasonably  liigh-speed,  large  memory  computer. 

Review  of  Problem  Formulation 

The  problem  of  a two-dimensional,  periodic  wave  propagating  in  water  of  uniform  depth 
has  been  discussed  in  Section  II  of  the  main  body  of  this  report.  If  the  water  is 
incompressible  and  the  motion  irrotational,  then  the  following  boundary  value  problem  can 
be  established  for  an  “arrested”  wave  system. 

Differential  Equation  (DE); 


Boundary 

Conditions 


= 0 (M) 

Bottom  Boundary  Condition  (BBC): 

w = 0,  2 = -h  (1-2) 

Kinematic  Free  Surface  Boundary  Condition  (KFSBC): 

15  “ z = n (X)  (1-3) 

Dynamic  Free  Surface  Boundary  Condition  (DFSBC): 
n + ^ (u  - O*  -t  = Q»  z “ n(x)  (14) 

^Motion  is  periodic  in  x with  spatial  periodicity  of  the  wavelength,  L.  (1-5) 


Equations  (I-l  through  1-5)  represent  the  common  formulation  for  all  of  the  cla—irjl 
nonlinear  water  wave  problems  in  which  it  is  assumed  that  the  wave  propagates  without 
change  of  form  and  a reference  coordinate  system  has  been  chosen  that  travels  with  the 
wave  form.  For  a specified  wave  height,  water  depth  and  wave  period,  the  goal  then  is  to 
determine  as  exactly  as  possible  a solution  to  the  formulation. 
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Stream-function  Solution 

The  general  form  of  the  Stream-function  solution  is: 


+ 


NN 

I X(n) 


n=l 


sinh 


cos 


a-6) 


The  water  displacement,  t],  is  determined  by  setting  z = t}  in  Equation  (1-6) 


T 

n « r 


NN 


E »n  - E J 1^  O'  * 

n*i  ^ 


n) 


cos 


^ x|  (1-7) 


where  is  the  (constant)  value  of  the  Stream-function  on  the  free  surface.  The  velocity 
components  are  defined  by: 


u - 


(1-8) 


(1-9) 


In  continuing  the  quest  to  determine  a solution  that  satisfies  Equations  (I-l)  to  (1-5)  as 
faithfully  as  possible,  it  is  noted  that  for  arbitrary  values  of:  and  the  X(n)’s,  the 

Stream-function  solution  exactly  satisfies  all  of  the  requirements  of  the  formulation  except 
the  DFSBC,  Equation  (14).  All  of  the  effort  can  therefore  be  directed  to  determining  these 
“free”  variables  such  that  they  represent  the  specified  wave  height  and  also  “best”  satisfy 
Equation  (14).  The  approach  employed  is  numerical  iteration,  in  which  a trial  solution  is 
regarded  as  available  and  at  each  step  of  the  iteration;  the  “free”  variables  are  modified  to 
improve  the  solution. 

As  a preliminary  step,  an  error  is  defined  in  the  one  remaining  unsatisfied  boundary 
condition, 

J 

E = i I (Q.  - 5)"  (MO) 

J j-i  3 


where  the  Q-'»  represent  equally  spaced  (in  B)  values  of  the  quantity  in  Equation  (14), 
and  Q represents  the  average  of  the  Qy’s.  If,  for  example,  J = 41,  and  the  free  variables 


100 


could  be  adjusted  so  tliat  E was  verj'  small,  then  the  associated  solution  would  provide  a 
good  fit  ft)  the  complete  formulation  at  these  41  points,  and  computations  have  shown  that 
the  fit  or  other  phase  angles  would  be  comparably  good.  The  problem  therefore  has  evolved 
into  one  of  minimizing  the  total  error  E.  The  procedure  ust^d  is  a least -squares  procedure, 
which  rctjuires  formally  that 


3E 

3L 


0 


(Ml) 


11 

3X  (n) 


0 


(M2) 


(The  parameter  is  not  determined  by  the  least -stjuares  jirocedure,  but  is  selected  such 
that  the  mean  water  level  is  not  changed  by  the  other  variables  selected.  This  will  be 
discussed  later.)  Examination  of  Equations  (1-1 1)  and  (M2)  further  will  indicate  that  the 
usual  least -s<{uarcs  procedure  is  not  applicable,  because  the  error  is  not  defined  as  a 
quadratic  function  of  the  unknowns.  This  problem  then  falls  in  the  category  of  a nonlinear 
least -squares  jiroblem. 

The  problem  was  linearized  as  follows.  Suppose  that  at  tlie  iteration,  a trial  solution 
is  available.  The  objective  is  to  select  chang<'s  in  the  unknowns  such  that  the  errors  will  be 
reduced.  If  this  were  a linear  least -squares  problem,  only  one  iteration  would  be  required. 
Expressing  the  quantity  Q in  terms  of  small  changes  in  the  unknowns  (to  be  determined  at 
the  k*^  iteration). 


where 

3Q  _ 3Q  3n  , ao  3u  . 3Q  3w  

ax(n)  3n  3X(n)  3u  3X(n)  3w  3X(n) 


3Q  _ ap  3n  ^ 30  3u  3Q  3w  3Q  3C 
^ “ 3n  3L  3u  3L  3w  3L  3C  3L 


(1-15) 


where  the  dQ/3Tj,  3Q/3u  are  obtained  from  Equation  (14)  and  the  3j7/3X(n),  3u/3X(n), 
etc.,  are  obtained  from  Equations,  (1-7),  (1-8),  etc. 
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Rewriting  the  least -squares  procedure  in  terms  of  the  unknowns:  AL  and  AX(n) 


= 

3 AL 


0 


(1-16) 


3E 

aAX(nT 


0,  n = 1 NN 


(1-17) 


Equations  (1-16)  and  (1-17)  represent  a set  of  NN  + 1 linear  simultaneous  equations  in 
terms  of  the  NN  + 1 unknowns.  After  each  iteration,  the  water  surface  is  recalculated,  by 
iteration,  from  Equation  (1-7)  and  V't,  >8  redetermined  such  that 


L 

n dx  = 0 

0 


(1-18) 


which  can  be  expressed  in  integral  form  as: 

^ (h  ri)|  cos  x]  dx  (1-19) 

' 0 

where,  in  the  computations,  a Simpson’s  rule  approximation  to  Equation  (1-19)  is  used. 

One  complete  iteration  comprises  a simultaneous  solution  for  AL  and  the  AK(n)’s  and  a 
redetermination  of  Successive  iterations  involve  exactly  the  same  procedure,  and  the 
iterations  can  be  terminated  when  successive  reductions  in  the  error  E are  small.  Numerical 
instabilities  can  occur,  especially  near  breaking  wave  conditions,  and  one  effective  procedure 
in  these  cases  is  to  apply  only  a fraction  of  the  AL  and  AX(n)’8  specified  by  the 
least-squares  solution. 

One  final  comment  should  be  directed  toward  the  problem  of  establishing  the  desired 
wave  height.  Although  it  is  possible  to  develop  more  sophisticated  procedures  which 
converge  on  the  wave  height,  the  procedure  followed  here  was  simply  to  conduct  successive 
runs  until  the  wave  hei^t  was  within  an  acceptable  limit  (1  percent)  of  the  desired  height. 


^ = I 

L 


X(n)  sinh 
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APPENDIX  n 


DEVELOPMENT  OF  COMBINED  SHOAUNG-REFRACTION  COEFnCIENTS 
Introduction 

This  appendix  describes  briefly  the  method  employed  to  calculate  the  combined 
shoaling— refraction  coefficients. 

Background 

TTie  shoaling— refraction  coefficients  developed  are  valid  for  a bathymetry  characterized 
by  straight  and  parallel  bottom  contours  and  for  a wave  system  which  suffers  no  energy 
losses.  The  two  principles  employed  are  Snell’s  Law  and  the  concept  that  there  is  no  energy 
flux  across  a wave  ray,  see  Figure  ll-l. 

Snell’s  Law  governs  refraction  and  relates  the  wave  propagation  speed,  C,  to  the  wave 
direction,  a , 


sin  Oi 

Cl 


Consti 


sin  02 

C;, 


(D-l) 


in  which  the  subscripts  pertain  to  any  arbitrary  depths. 

The  requirement  that  no  energy  is  propagated  across  wave  rays  may  be  written  as: 

cos  cos  aj^  = Constz  (^-2) 

in  which  Fj-g  represents  the  energy  flux  per  unit  width  in  the  direction  of  wave 
propagation  and  the  cos  a term  represents  the  width  between  adjacent  wave  rays. 
The  Fj>£  term  could  be  expressed  as  the  product  of  the  wave  energy  density,  TE,  and  the 
group  velocity,  Cg,  although  this  will  not  be  helpful  in  the  effort  here.  For  linear  wave 
theory,  it  is  possible  to  separate  the  refi'action  and  shoaling  effects  because  neither  the 
celerity,  C,  (governing  refraction)  nor  the  group  velocity,  Cq  (governing  shoaling)  depend 
on  wave  height.  For  our  case,  inspection  of  Equations  (Il-l)  and  (11-2)  will  show  that  the 
two  phenomena  are  coupled  throu^  the  dependency  of  C and  Cg  on  the  wave  height. 
Method 

The  method  employed  here  utilizes  the  dimensionless  energy  flux,  (Table  XI, 
Item  5)  and  the  dimensionless  wavelength,  L'  (Table  XI,  Item  1),  where 

F'  - !P- 

TE  -iFT 
T 

T I _ L 

^ TgTVJiT 
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Figure  IM.  Definition  sketch  for  shoaling-refraction  considerations 


Equation  (Il-I)  can  be  rewritten  in  terms  of  the  dimensionless  quantities  as: 


2tt  sin  ai 
gTi  L' 


27r  sin  a? 

gT7 


= Consti 


However  since  the  period  is  conserved,  i.e.,  Tj  = Tj 


(II-3) 


sin  ai 


sin  a. 


Consta 


(U4) 


The  energj  flux  relationship,  Fs^uation  (11-2)  can  be  expressed  as: 


21*1  * 

271 

[LoJ 

L’  cos  a = Constj 


or  recognizing  that  the  |)erit)d  is  cons<*r>'ed 

. 2 


LoJ 


L'  cos  a = Consti, 


ai-5) 


Equations  (II4)  and  (11-5)  describe  the  shoaling-refraction  process  in  terms  of  available 
dimensionless  parameters,  and  were  solved  as  described  in  the  following  paragraphs. 

Solution 

It  was  found  convenient  to  characterize  a particular  incoming  deepwater  wave  by  the 
direction,  a„,  and  deepwater  steepness,  /L^.  The  problem  is  to  determine  wave 
steepnesses  at  other  relative  depths  h/L„  such  that  Equations  (II4)  and  (II-5)  are  satisfied 
recalling  that  L and  both  def)end  on  h/L^  and  H/L^.  For  each  relative 
depth,  h/L^,  four  values  of  L'  and  F^g  arc  available  (for  H/Hg  = 0.25,  0.5,  0.75,  and 
1.0,  c.f.  Figure  23)  whereas  a continuous  distribution  is  required  for  the  purpose  here.  For 
each  relative  depth,  h/L^,  continuous  distributions  were  obtained  by  fitting  straight  lines 
between  the  four  available  points;  for  H/Hg  = 0,  it  was  assumed  that  the  simple  linear  wave 
theory  applied,  see  Figure  II-2  for  an  example  for  h/L^  = 0.02. 

For  given  /L^  and  the  constants  in  Equations  (II4)  and  (II-5)  are  defined.  The 
wave  steepness  H/L^  and  direction  □ at  any  relative  depth  are  determined  by  iteration  of 
the  two  following  equations. 


k-H 


= sin-»  [(L')*" 


(n-6) 


Lo 


k+1 


f(Ho/Lo)*  I'D  cos  Ool** 


cos  o* 


(n-7) 


in  which  the  superscript  k + 1 denotes  the  (k  + 1)*^  iteration  and  applies  to  the 
improved  estimates  of  a and  H/L^.  Once  these  estimates  are  known,  the  parameters  with 
the  k subscripts  on  the  right  hand  sides  of  Equations  (II-8)  and  (11-9)  are  calculated  and 
improved  estimates  of  a and  H/L^  are  determined,  etc.  The  procedure  was  initiated  in 
deep  water  and  the  wave  steepness  and  direction  calculated  at  the  remaining  nine  values  of 
relative  depth  advancing  shoreward  or  until  breaking  was  indicated.  At  each  relative  depth, 
the  iteration  converged  very  rapidly  with  three  or  four  iterations  usually  cufficient.  For  the 
first  iteration  at  a relative  depth,  the  initial  value  for  wave  steepness  was  taken  as  the  final 
value  fo'  iie  preceding  (greater)  relative  depth. 

The  shoabng-refraction  results  are  presented  in  graphical  form,  for  =0°,  10°,  20°, 
40°,  and  60°  in  Figures  25,  26,  27,  28,  and  29,  respectively.  A description  of  these  tables  is 
presented  in  Section  IV  and  two  examples  illustrating  their  appUcation  are  given  in 
Section  V. 
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♦ APPENDIX  III 

SAMPLE  SET  OF  WAVE  TABLES  FOR  CASE  4-D 
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